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PREFACE 


Interest  in  and  concern  about  the  problem  of  climatic  variation 
are  growing,  because  of  the  realization  that  earth's  climates  are 
undergoing  changes  that  can  exert  a profound  and  irreversible  Influence 
on  human  life.  It  is  now  considered  essential  that  we  explore  the 
mechanisms  of  climatic  variation  by  developing  numerical  climate  models, 
with  special  emphasis  on  the  climatic  effects  of  boundary  layer  (BL) 
phenomena. 

Parameterization  techniques  are  used  to  incorporate  BL  processes 
in  atmospheric  general  circulation  models  (GCMs) , which  serve  as  tools 
to  study  climatic  changes.  But  it  appears  that  owing  to  preoccupation 
with  the  treatment  of  various  other  physical  and  mathematical  features 
of  GCMs,  BL  processes  have  not  received  the  attention  they  deserve. 

It  is  now  realized  that  more  research  is  necessary  in  this  direction, 
because,  on  the  time  scales  relevant  to  climatic  changes,  BI  turbulent 
exchanges  (in  addition  to  other  physical  processes)  are  very  important 
in  governing  the  evolution  of  large-scale  processes.  We  therefore  need 
to  be  able  to  test  systematically  the  elements  of  our  understanding  of 
parameterization  hypotheses. 

This  report  surveys  existing  BL  formulations  that  furnish  the 
basis  for  parameterizing  the  BL  in  atmospheric  circulation  models.  Be- 
sides describing  BL  theories,  the  report  discusses  problems  encountered 
in  applying  the  conventional  theories  to  low  latitudes  and  to  the  atmos- 
phere over  oceans.  A separate  chapter  is  devoted  to  the  determination 
of  BL  height,  because  recent  research  has  indicated  that  this  is  the 
basic  parameter  of  boundary  layer  parameterization. 

This  report  is  a technical  contribution  to  the  broader  activities 
of  the  Rand  Dynamics  of  Climate  Program,  sponsored  by  the  Defense  Ad- 
vanced Research  Projects  Agency  and  directed  to  the  systematic  study  of 
( climate  variations.  Two  related  Rand  publications  by  the  present  author 
are: 
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• > Nwteriaal  Experiments  on  the  Computation  of  Ground  Surfaoe 

\ Temperature  in  ar.  Atmospheric  Circulation  Model , R-1511-ARPA, 
i May  1974. 

• Parameterization  of  the  Planetary  Boundary  Layer  in  Atmos- 
pheric General  Circulation  Models — A Review,  R-1654-ARPA, 
March  1975. 


SUMMARY 


It  is  generally  recognised  that  the  boundary  layer  (BL)  is  one  of 
the  major  energy  sources  and  momentum  sinks  in  the  atmosphere.  Re- 
search on  BL  dynamics  therefore  has  an  important  bearing  on  many  atmos- 
pheric problems.  More  emphasis  is  being  placed  on  the  study  of  the 
interaction  between  large-scale  atmospheric  dynamics  and  BL  processes. 
That  Interaction  is  particularly  important  to  researchers  attempting 
to  use  dynamic  methods  to  study  the  more  comprehensive  and  pressing 
problem  of  climatic  variation.  The  problem  involves  the  most  important 
question  of  all:  the  predictability  of  climatic  change,  whose  determi- 

nation requires  the  use  of  dypamic  models. 

Atmospheric  general  circulation  models  (GCMs)  are  being  used  as 
research  tools  to  test  hypotheses  and  to  predict  climates.  The  last 
few  years  have  seen  concerted  efforts  to  improve  both  simulation  model- 
ing and  the  long-range  prediction  of  climate  and  its  variation.  A 
significant  phase  of  that  improvement  has  resulted  from  the  effort  to 
understand,  in  addition  to  other  processes  such  as  convection  and  radi- 
ation, how  and  to  what  extent  BL  turbulent  exchanges  govern  the  evolu- 
tion of  large-scale  procenses.  In  fact,  one  of  the  most  formidable 
problems  confronting  atmospheric  scientists  is  the  whole  question  of 
determining  how  smaller-scale  conditions  (such  as  BL  turbulent  fluxes) 
are  excited  by  large-scale  conditions  and  in  turn  alter  them,  i.e.,  the 
question  of  " parameterization . " 

This  report  surveys  the  parameterization  "theories"  used  in  variors 
atmospheric  models.  There  being  no  comprehensive  fundamental  theory  of 
parameterization,  various  semi-empirical  theories  have  been  used.  Essen- 
tially, they  try  to  relate  horizontal  stress,  heat  flux,  and  moisture 
flux  at  the  earth's  surface  and  upward  through  the  BL  to  the  external 
(free  atmosphere  and  underlying  surface)  parameters. 

The  so-called  K-theory  has  been  widely  used  for  parameterizing  BL 
processes  in  atmospheric  numerical  motels.  The  assumptions  are  related 
to  the  form  of  "internal"  quantities,  such  as  eddy  coefficient  as  a 
function  of  bulk  parameters  (stability,  for  example),  and  the  relevant 
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equations  are  solved  usin^  appropriate  boundary  conditions.  Though 
based  on  more  or  less  arbitrary  assumptions,  this  theory  has  enabled 
us  to  gain  considerable  insight  into  the  behavior  of  BL  processes  in 
the  atmosphere. 

The  similarity  theory  is  based  on  the  hypothesis  that  the  turbu- 
lent regime  is  unambiguously  defined  by  the  values  of  the  parameter  UA, 
the  friction  velocity;  g/0,  the  buoyancy  parameter;  Q,  the  vortica! 

heat  flux;  and  f,  the  corlolis  parameter.  The  last  one  tmAes  it- coan-' 

■s 

lngless  to  apply  this  theory  near  and  at  the  equator.  Generally,  the 
results  of  treatment  of  the  BL  by  similarity  theory  show  satisfactory 

agreement  with  observations.  However,  sizable*  discrepancies  occur  be- 

\ 

tween  the  results  of  different  authors  concerning  the  valuer  of  "univer- 

\. 

sal"  constants  and  "universal”  functions  of  the  theory.  The seV discrep- 
ancies arise,  because  the  theory  has  been  developed  for  homogeneous 
barotropic  and  steady-state  conditions,  which  are  rarely  found  in  the 
real  atmosphere.  The  basic  theory  has  now  been  made  more  general  by 
taking  into  account  the  nonstationarity  and  influence  of  baroclinicity. 
Another  generalization  has  been  made  by  using,  in  place  of  conventional 
scale  height  UA/f,  a BL  height  h that  varies  in  time  and  space. 

On  the  basis  of  observations,  it  has  been  found  desirable  to  de- 
termine h through  a rate  equation  for  unstable  convection  BLb.  The 
use  of  a rate  equation  is  based  on  the  entrainment  hypothesis,  which 
is  applied  at  the  top  of  the  BL.  So  far,  however,  rate  equation?  have 
not  been  considered  appropriate  for  stable  BLs.  The  entrainment  hypoth- 
esis has  been  developed  basically  for  cloud-free  BLs,  though  there  ha  e 
been  attempts  to  incorporate  effects  of  radiative  fluxes  on  BLs  with 
clouds  within  them.  To  apply  the  entrainment  concept  in  practice,  it 
is  necessary  to  use  a closure  assumption  that  essentially  relates  the 
heat  flux  at  the  top  of  the  BL  to  that  at  the  underlying  surface. 

It  is  difficult  to  treat  the  BL  at  low  latitudes,  including  the 
equator,  because  none  of  the  existing  theories  are  applies ble  to  those 
regions  owing  to  the  presence  of  f,  the  coriolis  parameter,  in  almost 
all  basic  formulations.  There  is  no  alternative  but  to  use  the  exist- 
ing schemes  in  some  modified  form,  however,  since  no  schemes  have  be«  i 
designed  specifically  to  handle  the  low-latitude  BL  problems. 
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The  treatment  of  the  marine  BL,  on  the  basis  of  existing  theories, 
also  presents  a peculiar  problem  related  specif icslly  to  :ha  determina- 
tion of  "roughness"  of  the  ocean  surface.  The  r* facta  of  ocean  currents 
on  momentum  sources  also  are  to  be  considered.  Inasmuch  as  the  oceanic 
surface,  through  its  supply  of  heat  and  moisture,  strongly  affects  atmos- 
pheric processes,  concerted  research  on  parameterislng  the  marine  BL  is 
urgently  needed. 

Existing  theories  have  been  useful  for  regional  and  other  restricted 
studies,  but  are  likely  to  be  inadequate  for  global  models.  It  is  hoped 
that  this  survey  will  stimulate  further  application  of  semi-empirical 
theories  of  the  BL  to  global  models.  This  may  be  achieved  by  putting 
together  all  available  parameterization  schemes  in  a GCM  and  then  esti- 
mating the  advantages  and  deficiencies  of  different  approaches  through 
systematic  numerical  experiments. 
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Chapter  1 

INTRODUCTION  AMD  BACKGROUND 


SIGNIFICANCE  OF  THE  ATMOSPHERIC  BOUNDARY  LAYER 
IN  ATMOSPHERIC  CIRCULATIONS 

As  our  knowledge  about  the  larger-scale  dynamics  of  atmosphere 
grows,  more  emphasis  is  being  placed  on  the  study  of  atmospheric  bound- 
ary layer  (BL)  problems,  with  a consequent  emphasis  on  the  interaction 
between  these  two  scales  of  processes.  This  is  particularly  important 
as  the  more  comprehensive  problem  of  climate  is  studied  by  dynamic 
methods.  BL  researchers  must  review  their  understanding  in  the  con- 
text of  this  interaction,  while  people  studying  larger-scale  processes 
must  think  of  the  BL  as  an  Inherent  part  of  the  problem.  The  point 
deserves  emphasis  because  there  has  been  an  apparent  trend  to  treat  BL 
development  as  a closed  or  quasi-closed  problem,  and  use  research  find- 
ings only  to  provide  lower  boundary  conditions  for  the  free  atmosphere. 
Furthermore,  many  large-scale  dynamictsts  hold  the  view  that  a suffic- 
iently detailed  model  of  atmosphere  can  generate  its  "own"  BL,  and  that 
it  is  not  necessary  to  Include  the  BL  (as  a part  of  the  troposphere)  as 
an  input  to  the  free  atmosphere.  It  therefore  will  be  judicious  to  apply 
a brake  to  these  trends  and  study  the  role  of  the  BL  in  the  atmospheric 
general  circulation  in  a composite  and  comprehensive  way. 

It  is  generally  accepted  that  the  BL  is  one  of  the  major  energy 
sources  and  momentum  sinks  in  the  atmosphere.  This  layer  plays  a vital 
role  in  the  exchange  of  momentum,  heat,  and  moisture  between  the  earth's 
surface  and  the  free  atmosphere.  Computations  by  Kuug  (1969)  suggest 
that  almost  one-half  of  the  atmosphere's  kinetic  energy  is  lost  (dissi- 
pated) in  the  BL,  and  Wilkins  (1963)  estimates  that,  in  relation  to  the 
entire  troposphere,  90  percent  of  the  energy  dissipation  occurs  in  the 
lower  km.  The  atmosphere  also  receives  much  of  its  sensible  heat  and 
virtually  all  of  its  water  vapor  through  turbulent  processes  in  the  BL, 
which  eventually  influence  the  formation  of  weather  disturbances.  The 
role  of  the  BL  in  cumulus  convection  associated  with  fronts,  tropical 
storms,  and  cloud  clusters  ha-3  also  received  considerable  attention. 
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It  is  generally  believed  that  cimilus  convection  is  primarily  Induced 
by  the  frictional  veering  of  wind  io  the  BL.  the  so-called  CISK  mech- 
an  a (Conditional  Instability  of  Second  Kind),  which  is  also  induced 
by  frictional  convergence  in  the  BL,  has  been  regarded  as  necessary  for 
the  Intensification  of  tropical  storms.  Over  the  oaeane  the  top  of  the 
BL  has  been  observed  to  be  close  to  the  base  of  cumulus  clouds.  Oceanic 
cumuli  owe  their  existence  to  subcloud  layer  convergence,  believed  to 
be  formed  by  mechanical  turbulence  within  the  BL. 

Turbulent  motion  in  the  BL  is  a very  effective  means  for  the 
transport  of  momentum,  and  this  momentum  flux  from  the  atmosphere  onto 
the  earth  has  considerable  influence  on  the  evolution  of  weather  systems. 
And  since  one  major  goal  of  meteorological  research  is  to  improve  long- 
term (a  week  or  more)  weather  forecasts  and  then  use  the  new  techniques 
to  study  olimatia  changes , it  is  imperative  that  we  acquire  a better 
understanding  of  BL  dynamics,  including  energy  dissipation  and  all 
fluxes.  Currently,  long-term  forecasts  are  based  on  subjective  (syn- 
optic) methods  and  operational  numerical  models.  The  former  tend  to 
become  unreliable  after  2 to  3 days,  and  the  latter  leave  c.uch  to  be 
desired  after  5 days  or  more.  These  deficiencies  are  direct  y related 
to  the  relative  time  scales  of  energy  exchanges  within  the  troposphere 
vis-il-vis  turbulent  fluxes  through  the  BL.  For  example,  the  reaction 
time  of  the  atmosphere  to  turbulent  fluxes  of  heat  and  moisture  is  about 
3 days,  to  kinetic  energy  dissipation  or  to  latent  heating  about  1 day, 
and  to  radiative  fluxes  about  one  week.  In  view  of  the  significant 
influence  of  non-adiabatic  effects,  which  are  related  to  turbulent 
fluxes  of  heat  and  moisture  through  the  BL,  we  can  expect  significant 
improvements  in  long-range  forecasting  for  two  days  or  more  only  if  BL 
processes  are  well  understood. 

Even  though  the  most  important  energy  transactions  take  place  at 
the  earth's  surface,  we  cannot  simply  prescribe  boundary  conditions  of 
surface  characteristics  such  as  roughness,  temperature,  and  moisture, 
and  the ii  compute  the  input  of  energy  to  the  atmosphere.  We  canno-  be- 
cause the  atmosphere  itself  reacts  on  the  surface,  and  therefore  partly 
controls  the  boundary  conditions  and  the  energy  transfers  that  occur 
there.  By  virtue  of- the  interaction  between  small-scale  and  large- 
scale  atmospheric  processes,  the  study  of  turbulence  in  the  atmospheric 


BL  is  an  essential  factor  in  studying  the  physical  principles  of  long- 
range  weather  forecasting  as  well  as  the  theory  of  climatic  change. 

PURPOSES  OF  THE  REPORT 

The  purposes  of  this  report  are: 

• To  present  an  up-to-date  survey  of  the  theories  that 
have  been  used  to  parameterize  the  eddy  fluxes  of 
momentum,  heat,  and  moisture,  with  their  limitations, 

• To  discuss  the  determination  of  the  height  of  the  BL, 
considered  to  be  the  basic  parameter  of  BL  parameter- 
ization, and 

• To  discuss  the  specific  problems  encountered  in  the 
treatment  of  the  BL  in  low  latitudes  and  over  oceans. 

DEFINITIONS  AND  GENERAL  CHARACTERISTICS  OF  THE  BOUNDARY  LAYER 

The  BL  has  been  variously  defined  in  the  literature.  For  example, 
it  has  been  defined  in  general  terms  as: 

e That  portion  of  the  lower  atmosphere  in  which  the  wind 
deviates  from  gradient  or  geostrophic  flow  because  of 
the  retarding  influence  of  surface  friction. 

e The  region  adjacent  to  the  earth's  surface  where  small- 
scale  turbulence  is  induced  by  wind  shear  and/or  thermal 
convection  and  occurs  almost  continuously  in  space  and 
time  (Deardorff,  1972a).  (Small-scale  turbulence  is 
intermittent  beyond  the  BL.) 

e The  region  near  the  surface  in  which  turbulence,  of  a 
scale  not  much  greater  than  the  scale  height,  carries 
significant  fluxes  of  heat,  momentum,  water,  etc. 

(Charnock  and  Ellison,  1967). 

• A region  where  surface  effects  remain  important  but  are 
no  longer  completely  dominant  (Kraus,  1972).  This  defi- 
nition is  wider  than  generally  suggested;  it  includes 
regions  in  which  vertical  fluxes  of  momentum,  heat,  and 
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moisture  determine  the  vertical  distribution  of  these 
properties.  It  is  implied  thereby  that  there  are  as 
many  BLs  as  there  are  transported  properties.  However , 
the  various  fluxes  are  coupled  with  each  other  such 
that  the  BL  usually  can  be  represented  by  a single  layer 
in  which  surface  effects  remain  significant. 

As  these  definitions  indicate,  the  term  "boundary  layer"  applies 
to  a layer  of  air  above  the  earth's  surface  in  which  significant  fluxes 
of  momentum,  heat,  moisture,  and  matter  are  transported  by  turbulent 
motions.  This  definition  allows  for  the  inclusion  of  transport  brought 
about  by  penetrative  (cloudy)  convection,  which,  in  the  limit,  could 
apply  to  cumulonimbus  clouds.  Thus  turbulence  does  occur  in  cumulus 
clouds  and  along  frontal  surfaces  that  separate  air  masses.  But  since 
these  features  are  considered  to  be  mesoscale  explosions  (through  t..e 
BL)  into  the  free  atmosphere,  and  are  more  random  and  sparsely  scattered, 
the  turbulence  associated  with  them  is  usually  excluded  from  BL,  which 
(unlike  the  rest  of  the  atmosphere)  is  continuously  in  turbulent  motion. 

In  the  literature  the  most  commonly  used  terms  describing  the 
resolution  of  the  BL  into  various  layers  are 

• Interfacial  layer, 

• Surface  (or  constant  flux)  layer,  and 

• Ekman  layer . 

The  interfacial  layer  is  the  region  at  the  earth's  surface  that 
includes  the  immediate  neighborhood  of  a land  or  water  surface.  The 
processes  within  this  layer  are  not  well  understood,  particularly  at 
an  air-water  interface.  However,  considering  that  the  fluxes  through 
the  BL  must  be  estimated  from  conditions  in  the  interfacial  layer,  more 
research  is  urgently  needed  to  study  it.  For  example,  while  the  rough- 
ness parameter  over  land  has  been  assumed  to  be  independent  of  external 
parameters,  we  still  do  not  know  how  to  distinguish  between  aerodynamic 
roughness  and  topography.  This  knowledge,  if  available,  would  be  useful 
for  global  circulation  models,  which  have  to  consider  land  surfaces  of 


hit  ii  . 


-5- 


varying  complexities.  The  problem  of  specifying  the  roughness  param- 
eter for  ocean  surfaces  Is  still  more  complex.  In  this  case  the  rough- 
ness length  oarmot  be  directly  associated  with  any  geometric  parameter. 

It  has  been  suggested  that  it  may  be  related  to  the  slope  of  capillary 
or  short  capillary  gravity  waves  interacting  with  the  wind  (Kraus , 1972). 
Roughness  over  the  ocean  also  varies  with  the  fetch  and  direction  of 
the  wind. 

The  surface  (or  constant  flux)  layer  extends  from  immediately 
above  the  earth's  surface  to  a depth  of  some  tens  of  meters.  As  its 
name  implies,  the  vertical  fluxes  of  momentum,  heat,  and  moisture  are 
invariant  within  this  layer.  Though  based  on  inconclusive  observa- 
tional evidence,  "surface  layer"  has  also  been  defined  as  the  layer  in 
which  the  vertical  integral  of  the  time  variation  of  temperature  over 
the  layer  is  small  (<20  percent)  compared  with  the  magnitude  of  the 
heat  flux  (Lumley  and  Panofsky,  1964).  The  concept  of  constant  flux 
layer  (CFL)  is  useful  because  it  permits  us  to  estimate  the  value  of 
this  flux  from  observations  of  the  transported  quantity  at  varying 
distances  from  the  surface— for  example,  at  thi  surface  and  at  a level 
(usually  anemometer  level)  within  the  CFL.  In  fact,  this  layer  has 
been  studied  much  more  comprehensively  than  the  BL  as  a whole,  and  there 
now  exists  a reasonably  satisfactory  description  that  is  partly  theo- 
retical and  partly  observational.  The  transfer  properties  of  turbulence 
within  the  CFL  are  better  defined  for  near-neutral  and  unstable  condi- 
tions than  for  stable  conditions. 

The  Ekman  layer  is  thicker  than  the  constant  flux  layer  and  may 
fill  either  a part  of  or  the  entire  BL.  In  a classical  sense,  it  is 
the  region  where  the  vertical  flux  of  momentum  is  of  the  same  order  of 
magnitude  as  the  coriolis  and  the  pressure  gradient  forces;  this  essen- 
tially implies  a steady-state,  barotropie  situation.  In  the  real  atmos- 
phere, however,  the  conditions  are  often  infringed  by  synoptic  scale 
evolution  of  the  pressure  field,  as  well  as  by  the  diurnal  variation  of 
the  radiative  heat  flux,  both  of  which  alter  thermal  stratification 
and  thereby  affect  turbulent  fluxes  and  the  wind.  Nevertheless,  though 
the  classical  Ekman  BL  is  rarely  observed  in  nature,  its  concept  may 
be  applicable  for  certain  conditions  such  as  the  presence  of  a dense 
cloud  cover,  which  greatly  reduces  diurnal  variations  in  the  BL. 


BASIC  EQUATIONS  USED  IN  BL  NUMERICAL  MODELING 

The  ensuing  description  follows  that  given  by  Estoque  (1973) . 
Numerical  models  enable  us  to  examine  the  relative  l^ortance  of  ex- 
ternal factors  and  Internal  processes  in  determining  the  behavior  of 
the  planetary  BL.  The  external  factors  or  parameters  that  must  be 
considered  are  the  large-scale  synoptic  conditions  and  the  properties 
of  the  underlying  surface  of  the  earth.  The  internal  processes  are 
the  various  transport  processes,  the  most  important  being  the  eddy 
transport  along  the  vertical.  In  general,  the  large-scale  synoptic 
condition  is  specified  only  in  terms  of  the  geoatrophlc  wind,  G(z) 

(or  the  corresponding  large-scale  horizontal  pressure  gradient  Vp  ); 
in  some  cases  it  is  necessary  to  specify  the  large-scale  distributions 
of  potential  temperature,  &L(z),  and  the  mixing  ratio,  qL(z).  Terrain 
properties  considered  are  the  roughness  parameter  (zq)  , the  temperature 
(Tq),  and  the  mixing  ratio  (qQ) , all  of  which  may  vary  along  the 
horizontal. 

The  basic  problem  in  numerical  modeling  is  to  determine  the  space- 
time  variations  of  the  variables  (wind,  temperature,  and  moisture)  that 
describe  the  BL  as  functions  of  the  different  external  parameters.  So 
far,  generally,  the  two-dimensional  problem  (no  dependence  along  one 
horizontal  coordinate)  has  been  treated  more  extensively  than  the  three- 
dimensional  problem.  In  mathematical  terms,  a rather  general  formula- 
tion of  the  two-dimensional  problem  is  as  follows: 

GIVEN:  G(z)  or  VpL»  0L(z),  qL(r) 

z (x),  T (x,  t),  q (x,  t) 
o o o 

TO  FIND:  $(x,  z,  t),  w(x,  z,  t),  9(x,  z,  t) 

T(x,  z,  t),  p(x,  z,  t),  q(x,  z,  t) 

Here,  the  subscript  "LM  denotes  the  value  associated  with  the  large- 
scale  synoptic  flow  patterns,  while  the  subscript  "o”  refers  to  the 
value  at  the  earth  surface.  Each  unknown  variable,  which  may  be  re- 
garded as  the  8 urn  of  the  large-scale  value  plus  a perturbation  induced 
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by  BL  processes,  must  be  calculated  In  the  layer  between  the  surface 
and  a height  h.  Ideally,  the  value  of  h,  which  ve  define  to  be  the 
height  of  the  BL,  should  be  infinitely  large.  For  most  practical  pur- 
poses, however,  it  is  sufficient  to  assume  that  h is  of  the  order  of 
one  kilometer. 

The  set  of  two-dimensional  equations  that  must  be  solved  are: 


du  du  du  . , 1 dp.  1 ^Tx 

Jt  * u 3x  w 3z  v p dx  p dz 


dv  dv  dv 

dt  “ " u dx  “ w d* 


f„  _ 1.  & + 1 Hi 

ru  p 3y  + p dz 


(1.1) 


(1.2) 


30  30  31.3 

dt  " ~ u dx  “ w dz  + da  (V  + se 


(1.3) 


« — u — w + t — (F  ) + S 
3t  dx  dz  dz  ' qJ  q 


(1.4) 


du  . 3w  n 
dx  + * 0 


(1.5) 


lE  « E& 

3z  ■ RT 


(1.6) 


-&)1 


(1.7) 


These  equations  are  based  on  the  assumption  that  variations  of  dependent 
variables  along  y are  zero.  Most  of  the  notation  is  customary.  The 
quantities  t^,  t^,  Fq,  and  represent  the  eddy  fluxes  of  momentum, 
heat,  and  water  vapor;  Sq  and  are  the  internal  sources  of  heat  and 
water  vapor.  The  density  in  the  pressure  gradient  term  is  assumed  to 
be  constant.  If  the  eddy  fluxes  and  the  internal  sources  are  specified 
in  terms  of  the  unknown  variables,  then  the  equations  form  a closed  set. 

Boundary  conditions  are  to  be  specified  at  the  earth  surface  (as- 
sumed to  be  effectively  at  z * zq),  at  uPPer  boundary,  and  at  the 
lateral  boundaries.  The  usual  lower  boundary  conditions  are  as  follows: 
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Chapter  2 

APPLICATION  OF  K-THEORY  TO  BOUNDARY  LAYER  PARAMETERIZATION 

Atmospheric  perturbations  depend  on  spatial  and  temporal  varia- 
tions of  heat  and  moisture  flux  from  the  underlying  surface.  The 
relevant  transports  across  the  BL  are  carried  by  finite-amplitude, 
three-dimensional  perturbations.  Since  it  is  not  practicable  to  deal 
simultaneously  with  BL  perturbations  and  the  large-scale  circulations 
of  the  atmosphere,  the  former  effects  therefore  are  expressed  as  func- 
tions of  parameters , which  in  turn  can  be  determined  from  a knowledge 
of  large-scale  motions.  According  to  K-thecry,  this  can  be  achieved 
by  replacing  molecular  diffusivity  with  a stipulated  variable  eddy 
viscosity , K,  in  the  classical  laminar  flow  equations. 

For  parameterizing  the  BL  (as  per  K-theory  concepts) , the  atmos- 
phere is  usually  divided  into  two  layers:  (a)  a surface  (or  constant 

flux)  layer  and  (b)  the  so-called  Ekman  layer.  Parameterization  thus 
involve*  first  modeling  the  fluxes  at  the  surface,  and  then  the  turbu- 
lent diffusion  of  properties  (from  the  top  of  the  surface  layer)  through 
the  upper  Ekman  layer. 

TREATMENT  OF  THE  SURFACE  LAYER 

The  surface  layer  usually  extends  from  the  top  of  roughness  ele- 
ments (z  * z ) to  the  top  of  the  constant  flux  layer  (z  ■ z ) , which 
o s 

usually  represents  the  lowest  explicit  model  level.  The  quantity  z 

& 

can  take  on  values  up  to  a few  tens  of  meters;  it  la  usually  made  to 
coincide  with  anemometer  level  (~  10  m)  and  is  kept  fixed.  Estoque 
(1963)  used  a value  as  large  as  50  m for  z . However,  Sasamori  (1970) 

S 

used  variable  depths  for  the  constant  flux  layer  from  the  formula: 

<y2  • t„«>*  «-i> 

s 

where  T ia  the  response  time  and  K is  the  diffusion  coefficient  at 
M 

z * z ; the  formulation  shows  that  depth  is  large  when  the  intensity 
s 

of  transport  is  strong. 
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')•  fundamental  aquation*  witb  ’n  tha  layer  of  thicknasa  t are: 

tKfi)-) 

^(■SiH)-0/  <2-2> 

° ) 


whara  K^,  K^,  are , respectively,  eddy  diffusion  coefficients  of 
momentum,  heat,  and  moisture.  Busingar  at  al.  (1971),  Yamamoto  and 
Shimanukl  (1966),  and  Estoque  and  Bhumralkar  (1969)  have  given  formu- 
lationa  that  determine  tha  velocity,  temperature,  and  moisture  profiles 
in  the  surface  layer  over  a vide  range  of  diabatic  conditions.  These 
have  been  used  in  various  studies  of  the  atmospheric  BL,  both  local 
as  veil  as  general  circulation.  Notable  among  these  are  Sasamori  (1970) 
(microscale),  Delsol  et  al.  (1971)  (GCM) , Deardorff  (1972a)  (GCM) , and 
Pielke  (1974)  (mesoscale  sea  breeze).  Formulations  l>y  Businger  et  al. 
(19/1)  are  based  on  a large  number  of  direct  measurements  of  both  heat 
and  momentum  fluxes  under  both  stable  and  unstable  conditions.  Their 
formulations  are: 


and 


**  M _ a _ i5n~* 
U*  3*  u 


unstable  case 


(2.3a) 


— U*  — - — (1  - 95)"15 


«.**>  3‘ 


and 


U - (1  + 4.70  | Stable  case 


kz  u*  ae 


OW  9* 


(2.3b) 


U^0A  may  be  replaced 


by  (-v'9*) 


The  raclo  Kg/Ky  equals  0,74  and  is  baaed  on  observations  of  Bualnger 
at  al.  (1971).  In  the  above  expressions, 


<«  + « > 

£ where  L 


Deardorff  (1972a)  suggests  that  Eqs.  (2.3a)  and  (2.3b)  are  more  con- 
venient to  integrate  in  terns  of  familiar  functions  than  the  so-called 
KKYPS  (Luuley  and  Panofsky,  1964)  formulation. 

Yamamoto  and  Shimanuki  (1966)  developed  velocity  and  temperature 
profiles  in  the  surface  layer  under  various  stability  conditions.  The 
momentum  exchange  coefficient  near  the  ground  is  given  by 

kU„* 

*m " ♦jdci)  (2,4 

where  4^  is  a nondimensional  wind  shear  given  by 
. 3U  ( 

. 38  1 * 1 w 2 
^i  " U*  I in  all  cases 


4^  ~ unstable  conditions 
$2  ~ stable  conditions 


4>t  is  a function  of  nondimensional  height,  £,  given  by 


where 


m 


T5T 


^ «W 


6 = mean  6 for  the  layer 
0 * empirical  constant  - (15  ± 3)  . 


(2.6) 
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The  values  of  UA  and  9#  are  determined  from 

u*  * Ic^TlTn'  -”c±<Ico|)> 

*nd  k(e  - 6(*o)) 

6*  * <G±<|C|)  - Gjd^D) 

vhere  and  Gj  represent  profile  functions  obtained  as  a function  of 
|C|*  i.e.,  U(»  ^u ? + v^).  The  nondlmensional  height  |£|  ir»  given  by 


i • 1,2  in 
all  cases 


(2.7) 


The  values  of  and  for  particular  values  of  |£|  have  been 
tabulated  by  Yamamoto  and  Shlmanuki  (1966),  and  relevant  values  for 
and  can  be  obtained  from  interpolation  formulas  developed  by 
Shlmanuki  (1969).  The  relative  error,  using  these  approximate  rela- 
tionships, does  not  exceed  0.4  percent  for  unstable  and  0.6  percent 
for  stable  conditions. 

The  formulations  for  moisture  profile  in  the  surface  layer  are 
analogous  to  those  for  temperature.  The  diffusion  coefficients  for 
6 and  q in  surface  layer  can  be  obtained  from  by 


where  the  variable  0 (taken  as  1.35  by  Businger  et  al.  (1971))  may  be 
obtained  in  terms  of  dimensionless  wind  shear  from 

6 ■ ^ ^ j-g  | ^ for  neutral  and  unstable  cases 

and 


B - 1 


for  stable  cases  . 
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Estoque  and  Bhumralkar  (1969)  computed  the  wind,  temperature,  and 
moisture  profiles  In  the  surface  layer  from  the  following  equations: 


(^)  - * V. 


o k 
*o  + IT 


{ z + z \ 

1 ) 

j*n| 

' 

i 

o 

\ Zo  ) 

(z  + Z \ 

( a + acs) 

k ) 

In 

I o ' 

l Zo  } 

| (1  + a.S)  J 

(2.8) 


where  S Is  the  stability  parameter. 


a_.  Is  an  empirical  constant,  and  the  other  variables  have  usual  mean- 
ing. These  expressions  incorporate  the  effect  of  thermal  stratifica- 
tion somewhat  differently  than  does  the  so-called  KEYPS  formulation 
(Lumley  and  Panofsky,  1964)  and  appear  to  fit  empirical  data  better 
over  a wider  range  of  dlabatic  conditions. 

In  all  the  formulations  described  above,  UA,  0A,  qA,  the  integra- 
tion constants,  are  related  to  the  fluxes  through  the  surface  layer, 
which  are  given  by 


“U* 


*H  It  " a*6* 


h It  " u*q* 


And  since  the  surface  layer  is  characterized  by  constant  fluxes  in  the 

vertical,  these  are  also  fluxes  at  z * z , the  top  of  the  surface  layer. 

s 

In  all  the  preceding  formulations,  z over  land  surfaces  should 

o 

be  presented  realistically  as  a function  of  essentially  aerodynamic 
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Toughness.  The  values  of  20  to  70  cm  for  most  land  surfaces  have  been 
suggested  by  Fielter  et  al.  (1971).  However,  over  water , zq  cannot  be 
associated  in  any  direct  way  with  the  mean  surface  height  or  some  geo- 
metric parameter.  Charnock  (1955)  has  suggested  that  zQ  is  a function 
of  surface  stress  and  thus 


(2.9) 


m being  a proportionality  constant.  Clark  (1970)  used  this  formula  to 
obtain  zq  by  using  m - 0.032.  and  he  stipulated  that 


(z  ) . 2 0.0015  cm 

o water 


(2.10) 


One  may  use  a graph  compiled  by  Deacon  and  Webb  (1962)  to  determine  a 

functional  form  of  z over  water. 

o 


TREATMENT  OF  THE  EKMAN  LAYER 

Specific  detailed  formulation  of  K for  the  portion  of  the  BL  above 
the  surface  layer  has  been  suggested,  among  others,  by  Estoque  (1963), 
Yamamoto  and  Shlmanuki  (1966),  Deardorff  (1967).  Zilitlnkevich  et  al. 
(1967),  Estoque  and  Bhumralkar  (1969),  and  O'Brien  (1970). 

All  K formulations  involve  the  assumption  that  turbulent  transports 
are  proportional  to  the  gradient  of  the  transported  properties.  For 
example,  K for  momentum  for  a neutral  condition  is  given  by 

*2[(&)2  + (£)T 

where  l is  the  so-called  mixing  length,  and  u and  v are  horizontal  com- 
ponents of  wind.  In  the  following  discussions,  and  (which  are 
not  Independent  quantities)  are  considered  to  be  related  to  each  other 
through  the  following  relation  suggested  by  Businger  et  al.  (1971): 


1.35  . 
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However,  because  there  is  considerable  controversy  about  this  value, 
it  should  be  treated  as  tentative.  Various  authors  have  discussed 
the  vertical  distribution  of  K.  Whereas  well-known  and  accepted  re- 
sults are  available  for  specifying  K(z)  in  the  surface  layer,  the 
prescription  of  K in  the  Ekman  layer  is  not  quite  clear,  particularly 
for  numerical  models  in  which  a large  range  of  stability  is  expected. 
It  is  generally  assumed  that,  within  the  Ekman  layer,  K increases  to 
a certain  maximum  and  then  decreases  in  the  upper  portions.  Estoque 
(1963)  used  this  type  of  vertical  distribution  for  K in  a successful 
numerical  study  of  the  BL.  His  results  compared  very  favorably  with 
observations  of  the  Great  Plains  Field  Program.  At  the  top  of  the 
surface  layer  (50  m) , K was  estimated  from  the  results  of  surface- 
layer  formulations;  above  this  it  was  assumed  to  decrease  linearly 
and  reach  zero  at  the  top  of  model  atmosphere.  The  main  drawback  of 
this  profile  for  K (in  the  Ekman  layer)  is  that  it  fails  to  take  into 
account  any  variable  factors,  such  as  lapse  rate  and  wind  shear,  which 
can  affect  K.  Yamamoto  and  Shimanuki  (1966)  suggest  that  K (given  by 
Eq.  (2.4))  holds  throughout  the  Ekman  layer  and  thus 


9*Vl  8 


(S  - (?) 
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for  unstable  conditions 


for  stable  conditions 


(2.12) 


(See  Eq.  2.6.) 

Here  cr  is  an  empirical  constant  (Eq.  2.6)  and  p ■ 1/6  is  another  empir- 
ical constant  introduced  by  Yamamoto  and  Shimanuki.  Estoque  and 
Bhumralkar  (1969)  have  suggested  the  followir*s  generalized  version  of 
an  expression  devised  by  Blackadar  (1962)  •. 
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for  unstable  case 


for  stable  case 


(2.13) 
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where 


k(z  + z ) 
o 


k(z  + z ) 
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1 + 


kf (z  + zq) 
27  x Kf5  U 
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k(z  + z ) 
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(2.14) 


and  S,  the  stability  parameter. 
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a Is  an  empirical  constant  (see  Eq.  (2.8)).  Zilitinkevich  et  al. 
c 

(1967),  on  the  basis  of  dimensional  argument,  suggest  that 


where 


kz 


i+r’ 


A (the  asymptotic  value  of  2 at  large  z)  - 


(2.15) 


U 


JL 


f x 3.7  x 10' 


(2.16) 


It  may  be  noted  that  the  expressions  for  l in  Eqs.  (2.14)  and 

(2.16)  are  similar  except  for  some  slight  variations  in  the  constant. 
It  is  also  clear  that  the  effect  of  stratification  on  K (Eqs.  (2.13) 
and  (2.15))  arises  primarily  from  the  buoyancy  term  and  not  from  l, 
which  is  for  the  neutral  case. 

O'Brien  (1970)  proposes  a functional  form  of  the  eddy  coefficient 

K.  The  values  of  K and  their  derivatives,  with  respect  to  z (height 

8 

of  top  of  surface  layer)  and  h (height  of  BL) , are  used  to  derive  a 
Hermlta-interpolating  (third  order)  polynomial  given  by 
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K(z)  - K(h)  + 
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K(z  ) - K(h)  + 
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2(K(zs)  - K(h)) | 


h - z. 


(2.17) 


and  analogous  formulations  for  K for  6 and  q.  Here  it  is  assumed  that 


Wh  ’ 


0 i.e.,  the  eddy  flux  divergence  is  zero. 


The  values  of  and  at  zs>  the  top  of  the  surface  layer,  are 
determined  from 


kU.z 

* s 


w - $ (from  Eq*  (2*4)) 
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fcr  unstable  case 


for  stable  case. 


Here  nondimensional  height  | £ | is  given  by 

z 

8 


* 


where  L*  is  given  by  Eq.  (2.6). 

The  value  of  the  diffusion  coefficient  at  the  top  of  the  BL  (h) 
may  be  assumed  to  be  zero  (Sasamori,  1970)  or  be  arbitrarily  specified 
(Pielke,  1974).  This  is  done  essentially  to  prevent  significant  turbu- 
lent mixing  above  h.  Deardorff  (1971)  found  that  obtained  by  using 

Eq.  (2.17)  gave  too  large  values  at  intermediate  heights  compared  with 
the  values  in  the  surface  layer. 
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EFFECT  OF  THERMAL  WIND  AND  STABILITY  ON  WIND  PROFILE  IN  THE  BL 

Among  other  assumptions,  the  formulations  discussed  above  assume 
that  the  large-scale  pressure  gradient  is  constant  with  height. 

Blackadar  (1965) , Venkatesh  and  Csanady  (1974) , and  others  have  modi- 
fied the  barotropic  K-theory  to  include  the  more  commonly  prevailing 
baroclinic  conditions.  According  to  Blackadar  (1965),  to  incorporate 
thermal  wind  effects  it  usually  suffices  to  consider  pressure  gradient, 
which  varies  linearly  with  height.  In  terms  of  geostrophic  wind,  this 
Implies  that 

5 - K * 

where  £ is  the  surface  geostrophic  wind  and  Tu  is  the  vector  thermal 

O tl 

wind. 

Blackadar  (1965)  and  Blackadar  and  Ching  (1965)  have  obtained 
numerical  solutions  of  BL  equations  for  the  following  conditions: 

(1)  neutral  and  barotropic 
(ii)  stratified  and  barotropic,  and 

(Hi)  neutral  and  barcclinio. 

Figures  1 and  2 show  comparative  results  for  (1)  and  (ii) . One 
can  see  that  Instability  increases  the  surface  stress  and  decreases  the 
cross-isobar  inflow  angle  of  the  surface  wind.  Figures  3 and  4 show 
comparative  results  for  (i)  and  (ill) . The  differences  can  be  attrib- 
uted to  distortions  in  the  wind  spiral  (Fig.  5)  produced  by  the  thermal 
wind.  The  effects  of  stratification  are  more  pronounced  than  those  due 
to  thermal  wind,  e.g. , the  decrease  in  angle  a (from  neutral  condition) 
due  to  unstable  stratification  is  much  greater  than  that  produced  by 
thermal  wind  (parallel  to  surface  geostrophic  wind).  Also,  the  increase 
in  surface  stress  due  to  unstable  stratification  is  relatively  greater 
than  that  due  to  thermal  wind  alone. 

GENERAL  COMMENTS 

Questions  have  been  raised  about  the  treatment  of  the  atmospheric 
BL  as  consisting  of  the  surface  (constant  flux)  layer  and  the  Ekman 


Log  10  Ro 


Fig.  1— Graph  of  u*  versus  log  to  Ro  for  the  barotroplc  solution 
Effect  of  instability  Is  to  increase  the  surface  stress 
(Blackadar  and  Ching,  1965) 
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Log  to  Ro 


Fig.  2 — Graph  of  a versus  log  to  Ro  for  the  barotroplc  solution. 
Effect  of  instability  is  to  decrease  the  cross-isobar  angle 
(Blackadar  and  Ching,  1965) 
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for  thermal  wind  parallel  (dashed)  and  opposite  (dotted)  to  the  surface  geostrophic  wind 

(Blackodar,  1965) 
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A Leipzig  (Lettau) 

A O'Neill,  neutral 

□ Upavon  (Dobson) 

O Scilly  Is.  (Lettau) 

■I  Sheppard  and  Omar 

E = Eniwetok  Is. 
W = Wake  Is. 

J = Johnson  Is. 
SJ  = San  Juan 


X North  Sea  (Jeffries) 

• Braokhaven,  neutral  (Bernstein) 

0 Brookhaven,  baratropic  (Bernstein) 

■+■  Brookhaven,  baroclinic,  A =0,005 
(Bernstein) 

V U.S.  Great  Plains,  neutral  (Bernstein) 

• U.S.  vGreat  Plains,  all  cases,  barotropic 

O U.S.  Great  Plains,  all  cases,  A =0.005 
(Bernstein) 

▼ Haltiner  and  Martin,  neutral 


Fig. 4 — Angle  a between  the  surface  wind  direction  and  surface  geostrophic 
wind  direction.  Solid  curve  is  the  computed  relation  for  a barotropic 
neutral  layer.  Curves  are  also  shown  for  thermal  wind  parallel 
(dashed)  and  opposite  (dotted)  to  the  surface  geostrophic  wind 

(Blackadar,  1965) 
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layer.  Some  authors  have  suggested  that  the  BL  be  treated  as  a "single 
entity"  (Kraus,  1972;  Deardorff,  1973).  As  a consequence,  some  recent 
workers  have  developed  formulations  that  are  applied  to  the  entire  depth 
of  the  modeled  atmospheric  BL,  Instead  of  using  different  parameterisa- 
tions  for  different  regions  above  the  surface.  Orlanskl  et  al.  (1974) 
have  presented  a simple  formulation  for  K that  simulates  turbulent 
transfer  processes  over  different  layers.  Their  formulation  for  eddy 
coefficient  is  given  by 


K 


+ C 


(^>1 


1/3 


if  A6  < 0 (unstable  condition) 
II  A9  > 0 (stable  condition) 


Here,  K ,v  are  constant  values  of  eddy  diffusion  coefficient  and 
o o 

viscosity,  respectively,  and  Az,A8  are,  respectively,  local  values  of 
the  vertical  grid  size  and  potential  temperature  difference  across  the 
grid  box.  Thus  the  above  formulation  has  no  explicit  dependence  upon 
velociv.  'eformation  as  do  those  used  by  Estoque  (1963)  and  others. 

Here  it  is  assumed  that  turbulence  in  stratified  fluids  is  produced 
only  by  (i)  convective  instability  if  the  mean  stratification  is  un- 
stable, or  by  (ii)  wave  breaking  for  mean  stable  stratification.  In 
other  words,  it  is  implied  that  subgrid  turbulent  fluxes  occur  only  at 
those  locations  in  the  model  where  local  gravitational  instability  is 
produced  by  "resolved"  flow  processes  of  the  model.  Although  this 
formulation  does  not  involve  deformation  explicitly,  it  can  still  model 
turbulence  generation  caused  by  shear.  Also,  numerical  experiments  by 
Orlanskl  et  al.  (1974)  have  shown  that,  in  general,  there  is  no  con- 
stancy of  heat  flux  with  height  for  the  lower  50  m of  the  atmosphere, 
throughout  a daytime  simulation.  They  suggest  defining  "surface  layer" 
as  the  layer  in  which  vertical  integral  of  the  time  variation  of  temp- 
erature over  the  layer  is  small  compared  with  the  magnitude  of  heat 
flux.  They  did  find,  however,  that  a constant  flux  layer  can  be  justi- 
fied for  a few  hours  (2  to  3)  before  sunset.  They  also  found  that  the 
f height  of  the  constant  flux  layer  reaches  a minimum  after  sunrise  and 

| sunset  and  a sharp  maximum  (** 500  m)  late  in  the  afternoon. 


A 
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Yamamoto  et  al.  (1968)  have  also  sought  a formulation  applicable 
to  the  entire  BL.  As  described  by  Yamamoto  and  Shlmanukl  (1966) , the 
eddy  coefficient  K is  given  by 


K - 


4 


and  4 the  nondimens ional  wind  shear  by 


(2.18) 


♦4  - Id1'21*  43  - 242  + l - o 


where 


<*  f5-  <w 


It  may  be  noted  that  both  the  formulations  for  K and  K Involve  friction 
velocity  Ug  and  heat  flux  U*0*»  which  are  Independent  of  height.  Since 
this  situation  cannot  be  justified  observatlonally  in  the  upper  Ekman 
layer,  Yamamoto  et  al.  (1968)  have  suggested  that  if  U . is  replaced  by 
(x/p)  in  E q.  (2.18)  and  is  replaced  by  a prescribed  variable  heat 

flux  Q ■ "PCpK  36/8* , one  can  apply  the  resulting  formulations  to  the 
entire  BL,  as  in  Eqs.  (2.19)  and  (2.20): 
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(2.20) 


LIMITATIONS  OF  K-THEORY 

Though  still  widely  used,  various  problems  are  associated  with 
K- theory.  Some  of  these  are: 
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1.  Equations  (2.19)  and  (2.20)  Imply  essentially  that  the  sur- 
face layer  similarity  theory  has  been  extrapolated  to  cover  the  entire 
BL.  However,  since  shear  generation  of  Reynold's  stresses  can  become 
locally  insignificant,  it  is  not  always  advisable  to  modify  K-theory 
for  the  entire  BL.  For  example,  observations  have  indicated  that  shear 
generation  is  important  in  stably  stratified  low-level  flows,  but  there 
are  also  indications  (Lenchow,  1970)  that  it  can  vanish  in  a strongly 
heated  BL.  Consequently,  stresses  in  such  unstable  conditions  may  be 
almost  independent  of  small  local  mean  wind  shear.  In  such  cases,  K 
can  become  extremely  large  and/or  negative  aloft  in  the  BL  (Deardorff, 
1972b). 

2.  Another  important  drawback  of  K- theory  is  that  it  allows  only 
down-gradient  transports,  which  they  generally  are  not  when  flux  is 
effected  by  finite  amplitude,  relatively  large-scale  eddies.  Thus  the 
low-level  jet— the  velocity  maximum  within  the  BL — cannot  be  explained 
by  any  K-type  theory. 

3.  There  is  ample  evidence  that  under  conditions  of  upward  heat 
flux  the  lapse  rate  within  the  BL  is  slightly  less  than  adiabatic,  i.e., 
the  heat  flux  is  countergradient ; and  Lenchow  (1972)  has  demonstrated 
that  turbulent  momentum  fluxes  in  a heated  barodinic  layer  may  be 
locally  countergradient  above  a region  where  mean  wind  reaches  a maxi- 
mum. Deardorff  (1972c)  has  derived  a theoretical  expression  for  counter- 
potential temperature  gradient  that  can  sustain  an  upward  heat  flux. 

He  suggests  that  heat  flux  Q may  be  formulated  as 

<*  ■ -<Vh  (H  - *cg) 

rather  than  the  usual  form, 


* “ ‘‘Vh  fz 

Here  is  a small  positive  quantity  that  allows  a countergradient  up- 
ward  heat  flow.  His  findings  are  mainly  based  on  observations  of  Bunker 
(1956),  Telford  and  Warner  (1964),  Warner  (1971),  and  Lenschow  (1970). 
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Chapter  3 

SIMILARITY  THEORY  OF  THE  WHOLE  BOUNDARY  LAYER 


GENERAL  DISCUSSION 

When  K theories  were  found  rather  wanting  in  BL  studies,  atten- 
tion turned  to  determining  a direct  relationship  between  the  geo- 
strophic  wind  and  the  surface  stress  (or  friction  velocity,  UA)  without 
reference  to  any  hypothetical  eddy  coefficient.  These  efforts  led  to 
the  n ow  well-known  "similarity  theory,"  which  has  been  used  extensively 
In  BL  studies,  e.g. , Lettau  (1959),  Kasansky  and  Monin  (1960),  Zllltlnk- 
evich  et  al.  (1967),  Csanady  (1972).  Gill  (1968),  Blackadar  and  Tennekes 
(1968),  and  others. 

The  similarity  theory  essentially  recognises  that  the  details  of 
the  constant  stress  region  near  the  surface  cannot  be  modeled  by  the 
same  length  scale  that  characterizes  the  important  features  of  the 
planetary  BL.  In  other  words,  according  to  similarity  theory  the 
turbulent  BL  has  a distinct  double  structure  (Fig.  6),  consisting  of 
an  "inner"  (constant  stress)  layer  and  an  "outer"  (Ekman)  layer  (Clauser 
1956). 

The  similarity  theory  is  based  on  the  hypothesis  that,  for  hori- 
zontally homogeneous  and  stationary  conditions,  the  turbulent  regime 


Free  atmosphere 


h' 


Outer  layer 


* Zo, 


Inner  layer 


KIUMMflflKt  Surface 


Fig. 6 — Characteristics  of  the  boundary  layer 
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is  unambiguously  defined  by  the  values  of  the  following  "external" 
variables  (so  called  because  each  is  characteristic  of  the  location 
and  is  external  relative  to  the  regime  in  the  boundary  layer) : 

z “ rougnness  parameter  of  the  underlying  surface 
f * coriolis  parameter 
G * geostrophic  wind 
* buoyancy  parameter 

A0  = 6^  - 0q  ■ difference  of  potential  temperature  values  at  the 
upper  and  lower  limits  c,t  the  boundary  layer. 

The  internal  characteristics  of  the  bcttidary  layer  ares 

friction  velocity:  U* 

(or  shear  stress  at  surface) : tq 

angle  included  by  tq  and  Gq:  % 

the  surface  heat  flux:  Q (or  0*)  and 

o " 

stratification  parameter:  y. 

Each  of  the  external  variables  is  separately  a function  of  all 
the  internal  variables.  These  functional  relationships  are  best  formu- 
lated in  terms  of  the  following  nondimens ional  parameters: 


angle  a between  surface  stress  T and  G 
o o 


geostrophic  drag  coefficient  Cf 


PG2  'G  / 


geostrophic  heat  transfer  coefficient  • 


The  five  external  parameters  zq,  f,  G,  g/0,  and  A0  are  used  to 
form  two  nondimensional  groups  of  variables  such  as 
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Surface  Rossby  Number  = R ■ ■% — 

o fz 

o 


and 


Stability  parameter  S 


„ & M 

0 f G 


Empirical  evidence  (Deacon  and  Webb  1962)  indicates  that  water 
vapor  transfer  proceeds  by  much  the  same  physical  mechanism  as  heat 
transfer.  Thus  the  geostrophic  moisture  transfer  coefficient  may  be 
written  as 


E 


'Z 


pGAq 


where  E is  the  surface  moijture  flux,  and 
o 

Aq  » (q^  - qQ)  is  the  difference  of  mixing  ratio  at  the  upper 
and  lower  limits  of  the  BL. 


SIMILARITY  THEORY  AND  PARAMETRIC  RELATIONS 

According  to  similarity  theory,  the  dimensionless  velocity  and 
scalar  fields  in  the  BL  should  be  some  universal  functions  of 


a nondlmen8ional  height  £ 


•UL 


U, 


ku. 


and  a nondlmensional  stability  parameter  y = 


f L 


where  k ■ von  Karman  constant 


u: 


L ■ Obukhov  length  « - 


If]" 


The  above  similarity  parameters  are  essentially  chosen  on  the 
asstflp^lon  that,  under  all  conditions  of  stability,  BL  height  (h)  is 
proportional  to  the  height  scale  UA/f.  The  constant  of  proportion- 
ality is  considered  to  be  uniquely  determined  by  the  stability  param- 
eter y. 
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Let  us  scale  the  inner  layer  by  z and  the  outer  layer  by  U. /f . 

O " 

Let  us  also  assume  that  the  mean  velocity,  temperature,  and  moisture, 
together  with  their  vertical,  derivatives,  can  be  expressed  equally  well 
by  the  functions  appropriate  to  either  the  inner  or  outer  region.  Also 
assuming  horizontally  homogeneous  and  stationary  conditions,  and  after 
matching  the  similarity  profiles  in  a common  overlapping  region,  one 
can  get  the  following  relations: 
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(3.1) 


A8  1 [(  U* 


0, 


- C(y) 


- D(y) 


and  y components  of  the  surface  geostrophlc 

scaling  temperature, 
scaling  mixing  ratio, 

j is  the  surface  layer  value  for  the  ratio  of  the  eddy 
coefficients  for  heat  and  momentum  i>r  near-neutral 
conditions.  According  to  Buslnger  et  al.  (1971), 

Q|g  - 1.35  for  unstable  conditions. 

A(y),  3(y),  | are  the  universal  nondimens lonal  functions  depending 

C(y)  and  D(y) j on  the  internal  parameter  y - kUA/|f|L. 

Since  Eqs.  (3.1)  are  strictly  valid  for  barotropic  conditions, 
Ug0  and  V^o  may  be  replaced  by  and  V^,  the  geostrophlc  wind  com- 
ponents above  the  surface. 
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In  the  above  equations,  the  unknown  Internal  parameters  of  the 
BL  to  be  determined  are  the  friction  velocity,  UA;  angle  ag  between 
surface  shear  stress  and  geostrophic  wind  vectors;  6*,  the  temperature 
scale;  and  qA,  the  mixing  ratio  scale.  And  these  are  to  be  determined 
in  terms  of  known  external  parameters  f,  zQ,  buoyancy  parameter  g/6, 
Gq,  A6,  and  Aq.  The  free-atmosphere  (large-scale)  parameters  U,  V,  6, 
and  q are  considered  to  be  given. 

As  suggested  by  Monin  and  Zilitinkevich  (1967)  and  Clark  (1970), 
Eq.  (3.1)  can  be  inverted  to  obtain 


J?- vv  s>-k2{[tefe) 


a(S) 


-1 


+ B2(S) 


a " tan'1  u82  ■ W S) 

go 


In 


I# 


Sign  f 


A(S) 


(3.2) 


to . 


- C„(R  , S)  - 


G A6  H o 

O 


In 


w- 


C(S) 


E 

o 

GQAq 


W S> 


Equations  (3.2)  are  in  such  form  that  they  can  be  used  to  compute 
the  internal  characteristics  0*,  Q , and  if  the  external  character- 

istic8  of  the  BL,  zq,  f,  Gq,  A0,  and  Aq  are  given.  Here  CD»  CH,  Cg 
are,  respectively,  the  geostrophic  drag,  heat  transfer,  and  moisture 
transfer  coefficients.  Note  thu.t  in  Eq.  (3.2)  a bulk-type  stability 
parameter,  S - (g/0)/(A9/fG  ),  is  used  instead  of  y (of  Eq.  (3.1))  be- 
cause the  latter  contains  surface  fluxes  that  are  actually  being  sought. 
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MODIPICATIOM  OF  SIMILARITY  THEORY  EQUATIONS  (BPS.  (3.1)  AND  (3.2)) 

The  formulation*  described  above  (Eqs.  (3.1)  and  (3.2))  are  valid 
for  stationary,  horizontally  homogeneous,  diabetic,  and  barotroplc 
conditions.  These  enable  determination  of  surface  fluxes  in  terms  of 
the  surface  characteristics  and  the  mean  variables  evaluated  only  at 
the  top  of  the  BL.  Two  important  length- scales  are  involved:  the  BL 

height,  h (scaled  by  U^/f) , and  Monin-Obukhov  length  (L).  The  universal 
constants  A,  B,  C,  and  D are  dependent  upon  the  ratio  of  these  two 
length  scales.  Various  recent  studies  (Deardorff,  1972a,  1974;  Carson, 
1973;  and  others)  have  demonstrated  that  when  the  stability  parameter 
h/L  is  less  than  a v^lue  of  order  -1,  h is  not  related  to  U*/f.  As  a 
consequence,  Eqs.  (3.1)  and  (3.2)  have  been  modified  by  replacing  UA/f 
with  h,  the  actual  height  of  the  BL.  The  value  of  h itself  may  be 
determined  from  observations  or  from  a rate  equation  discussed  in  de- 
tail in  Chap.  4.  Arya  (1974)  has  independently  modified  these  equa- 
tions for  the  unstable  case  only,  because  of  great  uncertainty  in 
determining  h for  stable  conditions.  This  implies  that  the  parameters 
ought  to  be  switched  in  going  from  the  unstable  case  to  the  stable. 

One  can  avoid  that,  however,  by  using  h instead  of  U*/f  for  both  the 
unstable  and  stable  conditions  (Zilitinkevlch  and  Deardorff,  1974). 

Doing  so  is  not  likely  to  introduce  any  serious  errors  in  computation, 
particularly  considering  our  ignorance  of  how  h be  obtained  in  stable 
cases  (Helgarejo  and  Deardorff,  1974). 

In  view  of  the  above  discussion,  Eq.  (3.1)  can  be  written  by  re- 
placing U^/f  with  h,  the  height  of  the  BL.  Thus  we  have 
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where  Ui  ■ kh/L  and  subscript  i has  been  added  to  distinguish  universal 
constants  A,  B,  C,  D of  Eq.  (3.1)  from  those  in  Eq.  (3.3).  Note  that 
and  in  the  first  two  equations  of  (3.3)  are  reversed  in  defini- 
tion from  that  of  Zllitlnkevich  and  Deardorff  (1974). 

Similarly*  the  set  (3.2)  can  be  rewritten  (by  replacing  UA/f  with  h) 
as  follows: 
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Here  - ghA0/0G  is  the  stability  parameter.  Subscript  i distin- 
guishes the  constants  from  those  in  Eq.  (3.2).  may  also  be  obtained 
from  the  formula 
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It  may  be  recognized  here  that  the  stability  parameter  y - kU^/|f|L 
(of  Eq.  (3.1))  and  yt  - kh/L  (of  Eq.  (3.3))  are  not  Independent*  owing 
to  the  presence  of  L in  both.  Also,  under  relatively  stationary  condi- 
tions the  variation  of  L covers  a much  wider  range  than  that  of  U*/f 
or  h. 


Experience  confirms  Chat  Eqs.  (3.2)  and/or  (3.4)  are  more  con- 
venient to  use  in  practice  than  either  (3.1)  or  (3.3). 


EFFECT  OF  BAROCLINICITY  ON  WIND  PROFILES 
AND  GEOSTROPHIC  DRAG  LAW 

Equations  (3.1)  to  (3.4)  were  derived  for  barotropio  BL.  This 
is  perhaps  one  of  the  most  important  factors  that  cause  a large  scatter 
in  data  points  used  in  empirical  determinations  of  the  stability  de- 
pendent similarity  functions  A,  B,  etc.  (Aryat  1974;  Melgarejo  and 
Deardorff,  1974).  (The  other  factors  are  uncertainties  in  measure- 
ments and  effects  of  accelerations.)  Various  observations  have  shown 
a large  effect  of  barocllnicity  on  wind  shear  and  surface  cross-isobaric 
angle  aQ.  Some  theoretical  studies  also  have  studied  the  effect  of 
thermal  wind  on  wind  profiles  and  geostrophic  drag  laws  of  the  BL,  but 
most  are  based  on  K-theory  models  (Chap.  2),  which  suffer  from  the 
serious  limitation  that  they  do  not  consider  effect  of  geostrophic 
shear  on  values  of  K or  mixing  length.  Also,  K becomes  meaningless 
for  convective  conditions,  because  of  mean  gradients  being  close  to 
zero  or  of  the  wrong  sign  (Deardorff,  1972b  and  Wyngaard  et  al.,  1974). 

Following  similarity  arguments,  Hees  (1973)  modified  similarity 
functions  A ^ and  B^  (see  Eq.  (3.3))  to  incorporate  the  effect  of  thermal 
wind.  He  considered  the  special  case  of  geostrophic  shear  being  con- 
stant with  height  and  suggested  the  use  of 


in  place  of  A^(y^)  and  B^(y^),  respectively,  in  the  first  and  second 

equations  of  the  sets  (3.1)  to  (3.4).  Here  T and  T are  the  com- 

x y 

ponents  of  the  thermal  wind  invariant  with  height.  Arya  and  Wyngaard 
(1974),  through  a simple  physical  model  of  convective  BL,  have  general- 
ized the  above  concept  and  expressed  the  similarity  functions  Aio  and 
B^o  as  sums  of  a barotroplc  part  (A^  and  B^) , dependent  only  on  the 
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stability  and  BL  height  parameters,  and  a baroclinicity  dependent  part 
(A^(  B^),  by  using  the  "more  convenient"  surface  geostrophic  wind.  Thus 
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from  the  first  two  equations  of  (3.3).  According  to  Arya  and  Wyngaard 
(1974)  Ap  due  to  thermal  winds  are  given  by 


A’  * aM  cos  8 
i o o 

B*  - bM  sin  8 
1 o o 

where  a and  b ar.  -.oeffici-  s that  depend  on  the  variation  of  geo- 
strophic shea,;  if-h  height,  and  8Q  is  the  angle  between  the  surface 
geostrophic  shear  and  the  surface  wind. 

The  results  of  Arya  and  Wyngaard  (1974)  show  that*. 


e Actual  wind  shears  in  the  bulk  of  the  mixed  (convective) 
layer  are  small  compared  with  the  Imposed  geostrophic 
shear.  It  is  incor,  . rt  to  ...  ermine  geostrophic  shear 
simply  from  the  measured  wind  profile  in  the  upper  part 
of  the  BL. 

s The  flow  in  the  barcclinic  nr'  layer  is  far  from  geo- 
strophic. The  deviations  f-  ..he  geostrophic  equilib- 
rium are  most  striking  just  below  the  inversion  base, 
and  any  approach  to  the  former  must  occur  in  the  stable 
layer  above. 
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• aQ,  the  surface  cross- isobar ic  angle,  increases  toward 
the  equator,  a finding  that  agrees  with  observations  but 
is  opposite  to  the  trend  predicted  by  the  drag  laws  of 
Zllitinkevich  et  al.  (1967) . This  variance  can  be  attri- 
buted to  the  use  of  JA/f  as  BL  height  scale  instead  of  h, 
the  height  of  the  lowest  inversion  base, 
s Since  the  model  does  not  consider  an  evolutionary  BL, 
the  entrainment  effects  at  the  top  of  mixed  layer  have 
been  omitted.  Consequently,  the  level  where  stress 
vanishes  and  geostrophic  flow  is  attained  must  be  larger 
than  the  estimate  for  h used  in  this  study.  If  layer- 
averaged  geostrophic  winds  are  used,  however  (as  they 
are  in  general  circulation  models) , then  use  of  and 
B^  is  considered  adequate  for  describing  both  barotropic 
and  barodlnic  cases. 

DETERMINATION  OF  THE  UNIVERSAL  SIMILARITY  FUNCTIONS  A.  B,  C. 

AND  D.  AS  FUNCTIONS  OF  STABILITY  PARAMETER  U 

The  universal  functions  are  determined  empirically  from  observa- 
tions, or  theoretically. 

Evaluations  of  the  universal  functions  A,  B,  C and  D over  wide 
range  of  statistics  have  been  provided  by  Monin  and  Zllitinkevich 
(1967) , and  Zllitinkevich  and  Chalikov  (1968) , who  used  the  Great  Plain 
data  (Lettau  and  Davidson,  1957).  Clark  (1970)  has  also  determined 
these  functions  on  the  basis  of  observations  from  south  Australia. 
Examination  of  these  evaluations  has  shown  that  there  is  too  much 
scatter  and  too  many  differences  to  permit  reliable  estimates  of  the 
so-called  universal  functions.  For  example,  Carson  (1972)  has  compared 
the  plots  of  C(y)  obtained  by  Monin  and  Zllitinkevich  (1967)  and  by 
Clark  (1970)  over  a common  range  of  y (Fig.  7).  While  Clark's  curve 
is  well  defined  for  y < -50  it  shows  considerable  scatter  for  y > 50. 

On  the  other  hand,  the  points  used  to  determine  "Russian"  C(y)  (y 
ranges  from  -50  to  +70)  show  a wide  scatter,  thereby  emphasizing  marked 
differences  between  the  two  curves  throughout  the  common  tinge  of  y 
that  is  important  for  estimating  heat  fluxes.  Thus  there  is  an  apparent 
anomalous  behavior  of  C(y)  for  relatively  small  values  of  y;  it  can  be 
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FIfl»  7 — Variation  of  similarity  parameter  C with  stability 
as  suggested  by  Monin  et  al.  (1967)  and  Clark  (1970) 
(Carson,  1972) 


attributed  to  the  violation  of  restrictive  conditions  of  similarity 
theory  such  as  site  inhouogeneities,  thermal  wind,  etc.  It  is  also 
possible  that  some  spurious  effects  may  be  caused  by  the  differences 
in  determination  of  the  surface  fluxes  and  external  variables,  espec- 
ially the  choices  of  mean  flow  variables  such  as  8^  at  the  top  of  the 
BL. 

The  functions  A(y),  etc.  of  Eq.  (3.1)  or  (3.3)  can  be  determined 

if  we  know  U.,  9*,  z , U and  V , 69  and  6q.  Whereas  the  empirical 
" " o go  go 

determination  of  A(y),  B(y),  A^l^).  and  B^Uj)  is  relatively  straight- 
forward, that  of  C(y)/Ci(Pi)  is  quite  sensitive  to  the  specification 

of  69  (•  9,.  - 9 ).  Here  9 and  9.  are  the  potential  temperatures  at 
ho  on 

z ■■  z and  z » h,  the  top  of  the  BL  (the  same  is  true  for  D(y)  or 
o 

DjChj))*  Zilitinkevich  and  Chalikov  (1968)  used  a fixed  height  of  1 km 
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for  the  top  of  the  BL  (Irrespective  of  dynamic  and  stability  consider- 
ations), and  Clark  (1970)  assumed  a level  at  the  maximum  In  the 
U-proflle  as  the  top  of  the  BL.  It  la  evident  that  the  former  tech- 
nique is  likely  to  overestimate  0h  (especially  if  the  Inversion  base 
lies  below  the  fixed  BL  top  of  1 km) , and  the  latter  may  underestimate 
6^  because  the  BL  depth  in  this  case  is  about  half  of  the  actual  BL 
thickness  (Businger  and  Arya,  1974). 

Recently,  Arya  (1974) , after  reanalysing  the  data  from  the  two 
sites  used  by  previous  investigators,  has  suggested  better  estimates 
of  the  universal  functions.  Figure  8 shows  a comparison  of  empirical 
and  theoretical  determinations  of  C(y) . It  may  be  seen  that  whereas 
C(y)  varies  rapidly  for  the  stable  (y  > 0)  case,  it  shows  no  non- 
monotonic behavior,  as  is  implied  by  the  results  of  Clark  (1970) 

(Fig.  7).  However,  there  is  still  considerable  scatter  in  the  data 
as  y increases;  this  can  be  attributed  to  uncertainty  in  determlng  U^, 
6*»  Ugo,  Vgo,  etc.  under  very  stable  conditions.  For  the  unstable 
(y  < 0)  case  and  for  -y  <50,  C(y)  shows  a continuation  of  the  trend 
noticed  for  the  stable  case.  For  -y  as  50,  however,  the  function  shows 
no  trend  at  all.  This  implies  that  either  the  function  approaches  a 
constant  value  for  moderate  to  strong  instability,  or  that  the  rele- 
vant similarity  parameter  for  the  range  is  \l±  - h/L,  not  y ■ kU*/|f|L. 
But  as  shown  by  a comparison  of  C(y)  for  y < 0 and  C^y^)  (Fig.  9), 
at  least  in  this  case  the  behavior  of  C does  not  seem  to  differ  signifi 
cantly  when  y is  replaced  by  y^.  This  result  has  also  been  obtained  by 
Clark  and  Hess  (1973),  and  thus  this  study  has  not  established  the  use- 
fulness of  h as  scale  height  for  all  unstable  conditions.  The  very 
poor  comparison  of  Csanady's  (1967)  results  (Fig.  8)  with  empirical 
and  other  theoretical  studies  can  be  attributed  to  his  use  of  a simi- 
larity function  that  is  not  dependent  on  stability.  This  also  caused 
him  to  obtain  unrealistic  results  for  the  surface  cross-isobaric  angle, 
which  (in  his  study)  increased  in  unstable  conditions  and  decreased  in 
stable  conditions — a result  quite  opposite  to  observations. 

Melgarejo  and  Deardorff  (1974)  have  determined  the  similarity  func 
tions  A(y),  B(y),  C(y)  by  using  the  modified  BL  similarity  theory  (Eqs. 
(3.3)  and  (3.4))  in  which  they  have  used  the  observed  BL  height,  h,  as 
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Ch«  physical  length  scale  Instead  of  the  conventional  U^/f  for  both 
stable  and  unstable  cases.  Their  results,  as  well  as  the  scatter  of 
data  for  stability  functions,  are  comparable  to  those  obtained  by 
conventional  analyses.  There  Is  some  uncertainty,  however,  about  the 
use  of  h as  a relevant  parameter  for  stable  cases  because  of  the  poor 
validity  of  similarity  theory  on  the  stable  side.  In  fact,  Melgarejo 
and  Deardorff  (1974)  think  that  it  does  not  matter  whether  h or  UA/f  is 
used  for  stable  situations.  For  unstable  conditions,  however,  h is 
considered  store  relevant  than  U^/f  because  then  the  external  variables 
(A0,  Aq)  can  be  obtained  realistically  and  applied  toward  the  deter- 
mination of  surface  fluxes.  This  is  an  important  aspect,  specifically 
with  respect  to  general  circulation  models  of  coarse  vertical  resolu- 
tion. It  is  also  suggested  that  one  may  use  mean  actual  wind  at  h 
rather  than  geostrophic  wind,  because 

e The  results  are  then  applicable  closer  to  the  equator,  and 

e The  procedure  for  obtaining  universal  functions  A and  B 
is  more  analogous  to  that  for  determining  C (and  D) . 

DETERMINATION  OF  GEOSTROPHIC  DRAG  (C^) , HEAT  (C^) , 

AND  MOISTURE  (C~)  TRANSFER  COEFFICIENTS 

An  examination  of  Eqs.  (3.1)  and  (3.3)  shows  that  surface  fluxes 
are  Implicitly  Involved  in  the  definitions  of  the  similarity  parameters. 
As  a result,  these  equations  are  not. quite  suitable  for  computing  sur- 
face fluxes  from  the  known  values  of  external  variables.  However, 

Clark  (1970)  has  shown  that  once  A(y),  B(y),  C(y),  etc.  have  been  eval- 
uated empirically,  then  an  inversion  of  Eqs.  ('>.1)  and  (3.3)  results 
in  relationships  (3.2)  and  (3.4),  respectively.  And  the  latter  are  in 
a form  suitable  for  determining  internal  BL  characteristics  UA,  Qq,  E_ 
from  known  external  BL  characteristics,  namely  zq,  f,  G,  A9,  Aq.  For 
this  purpose,  we  need  external  stability  parameter  S (or  S^) , and  sur- 
face Rossby  number  R (or  h/z  ) . Clark  (1970)  has  obtained  values  of 

o o 

Ua/G,  aQ,  8^/A9,  q*/Aq  as  functions  of  S and  Rq  and  plotted  them  in 
the  form  of  nomograrps.  By  interpolation  on  these  nomograms,  U.  and  a , 
and  hence  T , Q and  E , may  be  readily  computed.  Figure  10  is  an 
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Fig.  10  — Nomogram  of  as  a function  of  R0  and  S for  near-neutral 
and  stable  conditions  based  on  Wangara  data 
(Arya,  1974) 


example  of  a nomogram  for  the  near-neutral  and  stable  case  that  can  be 
used  to  obtain  surface  stzess.  It  Is  seen  that  with  increasing  sta- 
bility U^/G  decreases,  and  its  dependence  on  the  surface  Rossby  number 
weakens.  Figure  11  shows  a corresponding  plot  of  U^/G  against  for 
unstable  (convective  conditions).  These  estimates  are  not  quite  re- 
liable, however,  because  of  a large  scatter  of  empirical  parameters. 


LIMITATIONS  OF  THE  SIMILARITY  THEORY 

Tne  similarity  theory  in  its  original  form  was  based  on  the  follow- 
ing simplifying  assumptions.  The  BL  was  considered  to  be: 
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Fig.  ll- W ..  Fig.  10  but  lot  0,/Go  « « •“*««  1'/Z' 

y . - . /-mnriitian*  (Arvo,  lr7A) 


• horizontally  homogeneous 

• barotropic 

• neutral 

• steady  state » and 

• with  no  radiative  heat  flux. 


It  la  evident  that,  though  the  reaulta  baaed  on  the  theory  In  thla 
lure  provided  reaeonable  Inalght  Into  the  determination  of  the  varla- 
blilty  ol  BL  characterlstlce , the  theory  cannot  repreaent  the  real 
atmosphere  conditions.  Consequently,  the  similarity  theory 
have  been  modified  by  adding  complexities  to  the  determining  P rame  ere. 
V„r  evamnle.  there  are  now  formulations  that  include  the  effects  of. 
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• horizontal  inhomogeneity 

• thermal  wind 

• diabatic  stratification. 

Introducing  these  modifications  violates  the  restrictive  conditions 
of  theory,  and  there  is  considerable  inconsistency  in  the  results  based 
on  the  similarity  theory.  This  is  especially  true  over  land,  where 
there  are  rapid  temporal  changes.  For  example,  for  a diurnal  cycle, 

A0 (-  0.  - 0q)  is  usually  large  during  early  morning  hours  (up  to  about 
0900  hr)  because  of  still-persisting  remnants  of  nocturnal  inversion. 

This  results  in  small  values  for  Q,  which  cause  anomalously  large 
values  of  C for  small  negative  values  of  stability  parameters  p (Clark, 
1972).  The  theory  is  also  of  doubtful  validity  in  the  regions  of  strong 
horizontal  discontinuities,  such  as  coastlines. 

From  the  point  of  view  of  "universality"  of  the  theory  applications, 
the  most  important  drawback  is  the  choice  of  U^/f  for  height-scale  for  all 
stability  conditions.  It  is  obvious  that  this  theory  (under  this  stip- 
ulation) has  no  effective  use  in  general  circulation  models  and  1b  out- 
right meaningless  at  the  equator.  As  discussed  earlier,  however, 
valiant  efforts  have  been  made  recently  to  modify  the  theory  to  render 
it  usable  for  more  general  purposes.  For  example,  while  the  scale 
height  may  be  considered  to  be  U*/f  for  stable  conditions  it  is  redefined 
to  be  h (the  height  of  the  base  of  inversion)  under  unstable  conditions. 

In  some  cases,  doing  so  perhaps  has  remedied  the  large  values  of  C re- 
ferred to  above.  Also,  the  scale  height  U^/f  has  been  made  a function  of 
stability  with  a view  to  studying  diurnal  variations  of  BL  character- 
istics (Zllltinkevich,  1972;  Carson,  1973).  There  are  also  suggestions 
to  use  actual  wind  at  h rather  than  geostrophlc  wind.  However,  since 
the  formulations  in  this  case  are  based  on  steady-state  similarity 
theory,  they  cannot  adequately  simulate  the  "observed"  evolutionary 
nature  of  the  BL. 

Under  unstable  conditions  and  especially  for  marine  conditions, 
the  similarity  theory  has  been  considered  to  be  of  questionable  valid- 
ity. The  reason  is  that  marine  conditions  are  usually  characterized  by 
inversions  (Kraus,  1972),  and  there  is  an  identical  mixed  layer  for  heat 


and  momentum.  Consequently!  the  specification  of  A0,  the  potential 
temperature  difference  between  the  bottom  and  top  of  the  mixed  layer, 
is  ambivalent  in  these  circumstances. 

The  formulations  of  similarity  theory  use  a quantity  zQ  implying 
a local  roughness  parameter — on  a mlcrometeorologlcal  scale — that 
varies  sharply  from  point  to  point  over  natural  surfaces  and  is  de- 
termined from  towers.  This  zQ  is  not,  in  general,  identical  with  a 
roughness  parameter  that  should  be  used  in  the  formula  for  the  large- 
scale  features.  This  has  led  Fiedler  and  Panofaky  (1972)  to  define 
an  effective  roughness  length  for  incorporating  surface  friction 
into  large-scale  models  of  the  atmosphere.  They  suggest  values  for 

z of  0.42  m (for  plains),  0.99  m (for  low  mountains),  and  1.42 
© 

(for  high  mounted*'**). 


m 
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Chapter  4 

DETERMINATION  OF  THE  HEIGHT  (h)  OF  THE  BOUNDARY  LAYER 


There  is  increasing  evidence  that  the  height  h of  the  BL  is  the 
basic  parameter  to  be  used  in  any  BL  parameterization.  This  chapter 
discusses  the  quantitative  determination  of  h using  both  diagnostic 
and  prognostic  techniques. 

DIAGNOSTIC  METHODS 

In  a very  detailed  survey , Hanna  (1969)  has  discussed  various 
diagnostic  methods  of  estimating  h.  He  has  indicated  that  for  many 
years  h was  estimated  on  the  basis  of  the  classical  theoretical  models 
of  Taylor  (1915),  which  predicted  that  wind  velocity  approaches  geo- 
strophic  velocity  in  an  asymptotic  manner  as  heights  increase.  This 
implied  an  estimation  oc  h by  means  of  arbitrary  criteria,  such  as  that 
requiring  h to  be  the  level  at  which  the  shear  of  the  wind  speed  first 
vanishes.  While  Blackadar  (1962),  Lettau  (1962),  and  others  have 
solved  equations  of  motion-yielding  asymptotic  wind  spirals  whose  thick- 
nesses are  functions  of  arbitrary  empirical  constants,  h has  also  been 
estimated  from  observed  soundings;  for  example,  Charnock  and  Ellison 
(1967)  used  radiosonde  ascents  to  determine  both  h and  other  BL  char- 
acteristics. The  survey  also  summarizes  the  various  methods  for  esti- 
mating h,  together  with  the  tests  of  their  results  against  the  1953 
O'Neill  BL  observations.  Table  1 lists  some  of  the  methods  he  con- 
sidered for  determining  h. 

More  recently,  Charney  (1969)  has  determined  h on  the  basis  of 
dynamics  of  the  turbulent  motion  in  the  BL.  He  postulates  that  h is 
determined  by  Ekman  layer  instability:  the  eddy  viscosity  increases 

until  the  eddy  Reynold  number  is  reduced  to  a critical  value.  Thus 


where  f is  coriolis  parameter  and  Ke  is  the  effective  eddy  viscosity 
determined  on  the  condition  that  Reynolds  number  is  equal  to  a critical 


SOURCE:  Hanna  (1969). 
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value  Rc>  which  has  been  found  from  experiments  to  be  100  for  neutral 
conditions . 

The  formulations  mentioned  above  are  by  and  large  applicable  only 
to  neucral  (or  near-neutral)  conditions.  Clark  (1970)  and  Carson  (1971) 
studied  observational  data  by  using  a stability  dependent  formulation 
for  h given  by 

u* 

h - C — F(y) 

where 

w ~ TfTT 

is  a stability  parameter.  C is  a constant  of  about  0.3.  Zilitinkevich 
(1972),  using  the  relationship  between  h and  "effective"  eddy  viscosity, 
Kg,  has  suggested  that  with  increasing  instability  h increases  as  | u i ^ 
(where  y is  a stability  parameter) ; and  with  increasing  stability  h 
decreases  as  |y|  . Using  different  reasoning,  Businger  and  Arya  (1974) 
have  also  shown  that  for  stable  conditions,  h varies  as  |y|-J*.  Fig- 
ure 12  shows  h defined  in  various  ways  as  a function  of  stability. 

PROGNOSTIC  METHODS 

In  the  preceding  section  we  considered  the  purely  diagnostic  de- 
termination of  h.  It  has  been  found  that,  for  steady-state  and  neutral 
or  near-neutral  cases  (Hanna,  1969) , a large  scatter  is  obtained  when 
theoretically  derived  h is  compared  with  profile-estimated  depths  of 
the  BL  (Fig.  13).  Even  if  h is  made  stability  dependent,  the  formula- 
tions (Carson,  1971;  Zilitinkevich,  1972)  are  totally  inadequate  for 
representing  the  BL's  evolving  nature  in  real  atmosphere.  Over  land, 
daytime  observations  typically  have  found  h to  increase  from  a height 
on  the  order  of  100  m shortly  after  sunrise  to  a height  of  1 to  3 km 
by  late  afternoon.  If  < 0 (where  is  the  large-scale  vertical 
velocity  at  the  top  of  boundary  layer) , then  the  maximum  height  reached 
may  be  much  less.  Similarly,  over  land  at  night,  h typically  increases 
from  a height  on  the  order  of  100  m shortly  after  sunset  to  only  200  to 
500  m by  early  morning.  In  addition  to  observations,  some  recent 


Th«  boundary  layer  height  defined  in  various  ways  as  a function  of  stability 

(Busing*  and  Arya,  1974) 
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UQ*/f) x 10” 3 (m) 


Fig.  13 — The  observed  boundary  layer  thickness  Hj 
at  O'Neill  as  a function  of  the  length  U0*/r 
(Hanna,  1969) 


theoretical  studies  (e.g.,  Deardorff,  1974)  have  demonstrated  the  in- 
adequacy of  a diagnostic  formula  for  calculating  h.  These  studies 
have  also  pointed  out  that  over  land  h generally  seems  determinable 
from  a rate  equation  only.  Figure  14  shows  that  variation  of  h given 
by  a rate  equation  compares  favorably  with  the  observed  h,  whereas  the 
diagnostic  formula  is  not  at  all  adequate.  It  may  be  remarked  here 
that,  while  there  is  general  agreement  regarding  the  use  of  a rate 
equation  for  h for  unetabte/neutral  conditions,  there  is  some  doubt 
about  its  use  for  8table  conditions.  Thus,  though  Deardorff  (1971) 
suggests  that  the  height  of  the  nocturnal  (stable)  boundary  layer  should 
be  predicted  by  a rate  equation,  Zilitinkevich  and  Deardorff  (1974) 


Fig. 14 — Height  of  the  boundary  layer  ( Deardorff , 1974) 


believe  that  a rate  equation  for  the  stable  case  nay  not  be  appropriate 
We  consider  below  a rate  equation  for  h that  is  appropriate  only  for  a 
convectively  unstable  BL. 

Physical  Descriptions  of  Convectively  Unstable  Boundary  Layer 

The  strong  heating  of  ground  by  solar  radiation  causes  thorough 
convective  mixing  in  the  lover  part  of  the  atmosphere,  which  in  turn 
establishes  a mixed  layer  where  potential  temperature  is  virtually 
independent  of  height.  It  is  often  observed,  however,  that  the  upper 
region  of  the  mixed  layer  is  slightly  stable.  Next  to  the  ground  it- 
self a shallow  superadiabatic  layer  occurs  (s£  a few  tens  of  meters) , 
marked  by  large  vertical  shears  of  wind  and  temperature;  here  the  heat 
is  predominantly  transported  by  mechanically  Induced  turbulent  motions. 
Above  the  mixed  layer  a deep  nonturbulent  stable  layer  occurs.  Separat 
ing  the  two,  however,  is  a highly  undefined  interfacial  entrainment 
layer  that  results  from: 

(i)  Physical  overshooting  into  the  stable  layer  of  convective 
elements  that  originate  in  the  surface  layer  and  continu- 
ally Intrude  into  the  stable  layer,  and 
(ii)  The  entrainment  of  capping  stable  air  into  the  mixed  layer. 

Both  factors  cause  an  increase  in  the  depth  of  the  mixed  layer 
(l.e.,  of  h).  Figure  15  is  a schematic  representation  of  the  develop- 
ing convectively  unstable  BL.  The  interface  entrainment  layers  vary 
widely  in  depth  and  character.  They  are  usually  represented  by  step 
discontinuity  in  potential  temperature  A6  (in  the  6 profile)  at  t » h, 
the  top  of  the  convectively  unstable  BL. 

The  other  physical  processes  that  affect  the  growth  rate  of  h are: 

e Horizontal  advection  at  h 

e Large-scale  vertical  velocity  at  h 

e Radiative  heat  flux. 


-52 


Z 

+ 


Deep  non -turbulent 


Turbulent  free -convection 


ae 

dz 


=0 


layer 


/Surface  layer  -|^<0 

-/I 

^ Forced  and  mixed  convection 


-*•  6 


Fig.  15 — Schematic  representation  of  the  developing 
convectively  unstable  boundary  layer 
(Carson,  1973) 


Though  not  much  is  known  about  the  first  item,  which  is  perhaps 
more  predominant  over  oceans,  the  second  factor  (as  discussed  else- 
where) significantly  influences  BL  evolution.  A subsequent  section 
discusses  the  third  factor  (radiation),  which  is  especially  important 
when  the  BL  contains  fog  or  stratus  or  stratocumulus  clouds. 

Rate  Equation  for  h for  the  Unstable  Case 

Many  authors  have  attempted  theoretical  time-dependent  determina- 
tions of  h for  the  atmospheric  mixed  layer.  See,  for  example,  Ball 
(1960),  Lilly  (1968),  Lavoie  (1968),  Deardorff  (1972a,  1973a,  1974), 
Tennekes  (1973),  Carson  (1973),  Stull  (1973),  and  Randall  and  Arakawa 
(1974).  These  studies  have  pointed  out  the  impossibility  of  establish- 
ing a stationary  state  except  in  very  special  circumstances.  The  mixed 
layers  have  been  found  to  continue  to  deepen  unless  counteracted  by 
subsidence  or  by  very  special  distribution  of  heating.  If  an  equilibrium 
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is  established  temporarily,  its  characteristics  depend  on  previous 
history  as  well  as  instantaneous  conditions.  Therefore  h aust  be 
treated  as  time-dependent  in  the  unstable  case.  Lavoie  (1968)  sug- 
gests a prognostic  equation  for  h,  namely, 


He  assumes  that  the  BL  top  is  a material  surface,  thereby  implying 
that  there  is  no  entrapment  of  mass  through  the  BL  top  from  the  over- 
lying  atmosphere.  However,  as  shown  by  observations  (both  in  the  atmos- 
phere and  by  laboratory  experiments) , entrainment  is  critically 
important  in  determining  h,  and  thus  Eq.  (4.1)  has  only  limited  use. 

The  last  term  in  Eq.  (4.1)  was  added  to  prevent  the  development  of  a 
auperadlabatic  layer  at  the  top  of  the  BL  whenever  the  inversion  is 
wiped  out  by  heating  from  below.  It  is  necessary,  of  course,  to  ensure 
that  the  physical  process  represented  by  this  term  does  not  violate  the 
treatment  of  the  thermodynamic  energy  equation. 

We  now  consider  a rate  equation  for  h that  incorporates  the  ef- 
fects of  entrainment  but  not  a radiative  heat  flux.  The  equation  may 
be  written  as 


Sr  - wh  ♦ we  <4-2> 

where  is  the  large-scale  vertical  velocity  at  the  top  of  the  BL  and 
Wg  is  the  entrainment  velooity.  As  indicated  above,  entrainment  pro- 
vides an  important  physical  mechanism  that  leads  to  deepening  of  the 
BL.  The  determination  of  entrainment  has  been  treated  both  theoreti- 
cally (Lilly  (1968),  Betts  (1973),  Tennekes  (1973),  Leardorff  (1972a, 
1973a,  1974),  Stull  (1973)  and  Randall  and  Arakawa  (1974)),  and  observa- 
tionally  (Deardorff  et  al.  (1969),  Lenschow  1973)).  Table  2 shows 
formulations  used  in  various  theoretical  studies  of  the  height  of  the 
convectively  unstable  BL;  these  formulations  are  primarily  based  on 
the  so-called  "jump"  models.  It  is  seen  from  these  expressions  that 
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entrainment  rate  depends  upon: 

the  downward  heat  flux  at  the  top  of  the  BL 
(A8)^  the  Intensity  of  the  inversion 

Y+  lapse  rate  just  above  the  jump,  and 
h height  of  the  BL 


Each  is  discussed  below. 

Q^,  the  Downward  Heat  Flux  at  the  Top  of  the  BL.  This  is  an  im- 
mediate consequence  of  ‘-.he  entrainment  of  warm  air  into  BL.  The  physi- 
cal explanation  for  the  maintenance  of  downward  heat  flux  at  the  BL  top 
can  also  be  seen  from  a consideration  of  kinetic  energy  of  turbulence 
near  the  upper  limit  of  the  mixed  boundary  layer,  where  entrainment 
occurs.  In  general,  if  the  dissipation  rate  of  kinetic  energy  is  neg- 
ligible, there  must  be  a flux  convergence  of  kinetic  energy  to  maintain 
a downward  heat  flux.  Tennekes  (1973)  suggests  that,  on  the  basis  of 
dimensions  of  this  flux  divergence. 


g 


o 3 

a 

w 


(4.3) 


where  Ow  is  the  standard  deviation  of  vertical  velocity  in  the  BL  and 
g/0o  is  the  buoyancy  parameter.  Also,  for  a BL  in  a state  of  free 
convection,  when  the  turbulent  kinetic  energy  is  maintained  by  buoyancy 
only. 


(4.4) 


(Tennekes,  1970;  Deardorff,  1972a  and  others).  Combining  Eqs.  (4.3) 
and  (4.4)  yields: 


(4.5) 


a result  that  agrees  well  with  the  experimental  data  of  Deardorff  et  al. 
(1959)  for  A = 5.  Here,  parameter  A may  be  considered  the  measure 
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of  the  degree  of  interfacial  entrainment  at  h.  Equation  (4.5)  consti- 
tutes the  closure  hypothesis  for  the  entrainment  theory  at  the  top  of 
the  BL. 

Ball  (1960)  assumes  a value  of  unity  for  A which,  though  not 
realistic  in  the  light  of  the  experimental  results  of  Deardorff  et  al. 
(1969),  provides  a maximum  possible  entrainment  criterion  for  models 
driven  solely  by  surface  heating.  Lilly  (1968)  suggests  that  the  mag- 
nitude of  actual  downward  heat  flux  lies  somewhere  between  0 anJ  1. 

More  recently,  Deardorff  (1974)  and  Betts  (1973)  have  used  a value  of 
0.25  for  A in  their  models.  Carson  (1973a)  has  suggested  that  A should 
be  time-dependent — partly  determined  by  wind  shear  across  the  entrain- 
ment layer. 

(A8)^,  the  Intensity  of  the  Inversion.  This  represents  a positive 
jump  in  average  virtual  potential  temperature  at  h.  Lilly  (1968) 
suggests  that  the  entrainment  rate  is  inversely  proportional  to 
(A0)^.  For  a realistic  determination  of  W£,  however,  it  Is  extremely 
important  to  consider  the  time  variation  of  (A0)^  also.  Tennekes  (1973) 
uses  a rate  equation  for  (A0)^  to  account  for  kinematic  behavior  of 
(A6)^.  He  postulates  that  while  on  the  one  hand  (A0)^  tends  to  decrease 
as  the  BL  heats  up,  on  the  other  hand,  it  tends  to  increase  as  the  en- 
trainment proceeds  at  h.  The  resulting  equation  is 


_d 

dt 


(4.6) 


The  second  and  third  terms  combined  represent  the  situation  inside  the 
BL.  It  can  be  seen  that  this  use  of  Eq.  (4.6)  obviates  the  problem 
Inherent  in  using  the  formulation 


tf  ■ 

"e 


(4.7) 


to  determine  h,  as  (A0).  -*■  0. 

+ n 

Y , the  Lapse  Rate  of  Potential  Temperature  Jv " t Above  the  Jump 
at  h.  Deardorff  (1972a,  1974)  and  Tennekes  (1973)  have  suggested  that 


entrainment  rate  is  inversely  proportional  to  lapse  rate  in  the  stable 
air  above  inversion.  Deardorff  (1972a)  used  this  postulate  only  when 
(A0)^  approached  a certain  arbitrary  small  value.  In  a later  paper 
(1974),  he  regards  (A8)^  as  small  or  uncertain,  and  evaluates  entrain- 
ment by 


1.2  Q 


(4.8) 


Carson  (1973)  uses  y+  to  evaluate  the  intensity  of  inversion  as  a func- 
tion of  parameter  A (Eq.  (4.5))  and  of  h that  varies  with  time. 

Height  h of  the  Convectively  Unstable  BL.  If  h becomes  very  small, 
most  of  the  entrainment  formulations  listed  in  Table  2 break  down. 
Actually,  for  small  h,  the  dynamics  of  the  surface  layer  would  tend  to 
dominate  that  of  the  entire  BL.  As  a result,  it  has  been  considered 
advisable  to  prescribe  a lower  limit  for  h.  For  example,  Deardorff 
(1972a)  has  suggested 


hmin  ' za  + 50  z„ 

where  z is  the  surface  height,  and  z is  the  roughness  length.  Tennekes 
s o 

(1973)  puts  h . at  100  m. 
min 

"Interpolation"  Rate  Equations  for  Determining  Unstable 
Boundary  Layer  Height 

We  have  discussed  above  the  rate  equation  of  h for  a convectively 
unstable,  dry,  inversion-capped  BL.  The  physical  problems  associated 
with  entrainment  are  not  well  understood  as  yet,  and  hence  most  of  the 
treatments  are  exploratory.  Recently,  Tennekes  (1973)  and  Deardorff 
(1974)  have  suggested  rate  equations  for  h that  are  essentially  inter- 
polation formulations.  For  example,  consider  the  equation 
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Again,  if  both  and  (A0)^  tend  to  zero,  dh/dt  - becomes  Indeter- 
minate. However,  Deardorff  (1974)  found  from  his  numerical  experiments 

that  when  both  Q.  and  (A8).  become  zero, 
n n 


or 


f - wh  - 0.2  (f-  hQo) 


t ' Wh  ' °'2  M* 


1/3 


(4.10) 


where  the  quantity 


W*  S( 


8_ 

0 


hQ  i 
o °/ 


\l/3 


is  the  mixed  layer  convective  scale. 
He  derives  the  following  interpolation  formula  by  combining  Eqs.  (4.9) 
and  (4.10): 


1.8  Q 


dh 

dt 


*r-  W, 


1 + 1 


•'§K)1 


hy+  + 


(‘  * *'*  I) 


(4.11) 


Equation  (4.11)  prevents  the  development  of  a singularity  when 


(i) 

h ->  0 

or 

Y*  0 

(ii) 

Q 0 

o 

and 

Y+  0 

(iii) 

Q -*•  0 
o 

but 

yNo 

Deardorff  (1974)  tested  Eq.  (4.11)  (with  - 0)  against  3-D  model 
values  of  h as  well  as  against  observed  and  diagnostic  values  of  h 
(Fig.  14) . The  disagreement  between  observed  h and  that  obtained  by 
rate  Eq.  (4.11)  and  the  3-D  model  can  be  attributed  to  his  assumption 
of  ■ 0.  Empirical  adjustment  of  the  numerical  coefficients  in  Eq. 
(4.11)  may  be  necessary  to  obtain  accurate  results. 

Tennekes  (1973)  has  also  suggested  an  interpolation  formula  to 
treat  cases  in  which  BL  turbulence  is  maintained  by  buoyancy  as  well 
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as  mechanical  mixing.  Thus, 


dh 

dt 


e . 
2.5  -St  I 


u: 


8 h (A0). 


+ 0.2 


<A0>h 


(4.12) 


This  formula  cannot  be  used  if  h is  too  small.  Also,  for  (A8)h  "0, 
dh/dt  becomes  infinite,  a physically  untenable  situation.  In  that 
case  a limiting  value  given  by 

^ “ 0.2  (Ow)  (4.13) 

should  be  prescribed  for  the  entrainment  velocity.  Equation  (4.13), 
which  is  similar  to  Eq.  (4.10)  used  by  Deardorff  (1974),  gives  the  rate 
at  which  a BL  with  (A0)^  and  y ■ 0 entrains  aloft  (Tennekes  and  Lumley, 
1972). 


Rate  Equation  for  a Stable  Boundary  Layer  Height 

The  BL  is  almost  always  stably  stratified  over  land  at  night,  when 
it  is  typically  much  shallower  than  it  was  the  previous  afternoon.  As 
discussed  at  the  beginning  of  this  chapter,  h has  evolution  character- 
istics for  the  stable  case  also,  though  the  variations  are  not  as  pro- 
nounced as  convective  conditions.  It  was  also  pointed  out  that  there 
is  some  doubt  about  the  validity  of  using  a rate  equation  for  h under 
stable  conditions  because  BL  depth  then  is  usually  small  (in  comparison 
with  the  similarity  scale  height  0.3  U^/f)  and  not  likely  to  be  much 
affected  by  vertical  velocity  and  turbulence,  which  are  themselves  small. 
Deardorff  (1972a)  suggests  an  interpolation  formula  between  the  neutral 
height  of  BL  and  a value  of  h proportional  to  Monin-Obukhov  length  L. 

His  formula  is 

where  H^,  is  the  height  of  tropopause.  The  last  term  ensures  that,  for 
a neutral  condition  a*.  the  equator,  h does  not  exceed  H^. 


-60- 


Alternatively,  h can  be  set  at  50  m at  the  start  and  then  allowed 
to  Increase  at  a rate  proportional  to  UA  whenever  surface  flux  (or 
stability)  changes  sign  from  positive  .to  negative. 

Despite  an  uncertainty  regarding  use  of  a rate  equation  for  a 
stable  BL,  Deardorff  (1971)  did  determine  h for  a stable  BL  by  means 
of  a prognostic  equation: 


- Wh  - 0.025  U*[l  - h/ (0. 35  Vjf ) ] 


(4.14) 


The  factor  in  brackets  has  been  included  on  the  supposition  that  the 
similarity  theory  for  a neutral  case  provides  an  upper  limit  to  h in 
s table  atmosphere.  Thus 


§ - WH  - 0.025  U*  Max 


1 - 


V 

0.35  j- 


0 


(4.15) 


The  factor  (0.025)  in  this  equation  is  only  tentative  and  should  be 
revised  on  the  basis  of  observations.  For  example,  observations  of 
Clark  et  al.  (1971)  suggest  a value  of  0.04  or  0.05.  Randall  and 
Arakava  (1974)  have  used  another  version  of  Eq.  (4.15)  to  predict 
boundary  layer  height  for  stable  conditions: 


1 

Wh“  O'025  U*  ) 1 ‘ Ri”’  0 

c 


(4.16) 


where  Rig  is  bulk  Richardson  number  and  Rig  is  a critical  value  of  Rig. 
According  to  Deardorff  (1972a), 


Ric  - 3.05 


RADIATIVE  FLUX  EFFECTS 

Previous  sections  have  discussed  a rate  equation  for  h that  in- 
corporates effects  of  upward  eddy  heat  flux,  large-scale  vertical 
velocity,  and  entrainment  at  the  interface  between  the  mixed  layer  and 
the  overlying  stable  atmosphere.  These  discussions  did  not  consider 
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the  effect  of  radiative  heat  flux  on  the  evolution  of  h in  the  presence 
of  clouds,  though  the  presence  of  water  vapor  can  be  accounted  fo;  by 
using  virtual  potential  temperature. 

Lilly  (1968)  first  put  forward  a theory  that  considered  radiation 
off  the  cloud  tops  as  an  essential  element.  His  theory  was  applied  to 
a marine  boundary  layer  containing  fog  or  stratus  or  stratocumulus 
clouds  typically  observed  near  the  California  coast,  northern  Chile, 
southern  Peru,  and  southwest  Africa  most  of  the  year.  On  the  basis 
of  his  theory  he  concluded  that  the  observed  Inversion  (in  the  trade- 
wind  regime)  at  a height  of  500  to  1000  m and  with  an  intensity  of 
15  to  20°C  cannot  be  sustained  without  a radiatively  effective  cloud 
cover.  This  pointed  out  the  importance  of  incorporating  radiative 
heat  flu:i  effects  in  studying  the  evolution  of  a cloud-topped  BL.  Phys- 
ically, for  such  a BL,  strong  flux  divergence  gives  rise  to  strong 
radiative  cooling  (destabilization)  at  the  BL  top,  which  in  turn  causes 
a rise  in  BL  height.  Following  Lilly  (1968),  Randal]  and  Arakawa  (1974) 
have  described  the  buoyant  instability  of  a cloud-topped  BL  in  some  de- 
tail. They  have  derived  an  expression  for  entrainment  rate  that  in- 
corporates effects  of  radiative  heat  flux  as  well  as  the  presence  of 
liquid  water  within  the  BL.  However,  according  to  Lilly  (1968),  the 
choice  of  entrainment  approximations  does  not  appear  to  strongly  affect 
the  general  character  of  steady-state  solutions. 

Whereas  the  temporal  evolution  of  a clear  BL  is  fairly  well  known, 
that  of  a cloud-topped  BL  is  still  in  the  speculative  stage.  In  a 
cloud-free  BL,  for  example,  net  buoyancy  through  the  BL  and  kinematic 
surface  heat  flux  generally  have  the  same  direction  (Fig.  16) . The 
picture  is  less  certain  for  a cloud-topped  BL.  It  is  possible  that  in 
this  case  net  buoyancy  and  surface  heat  flux  may  act  oppositely  to  each 
other  (Fig.  17).  Net  buoyancy  of  the  BL  can  be  positive  (upward),  even 
though  surface  flux  is  negative  (u: inward)  (Randall  and  Arakawa,  1974). 
Lilly  (1968)  had  suggested  this  argument  on  the  consideration  that  while 
the  cloud  layer  is  maintained  in  an  active  turbulent  state  because  of 
the  release  of  latent  heat  there,  the  sub-cloud  layer  is  characterized 
by  a downward  heat  flux. 


Stratus  layer 


6 12  18  24 


N«t  buoyancy 


Surface  heat  flux  Qq 


Fig.  17  — Diurnal  cycle  of  cloud  topped  boundary  layer 
(Randall  and  Arakawa , 1974) 
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Mak  (1974)  has  also  studied  the  role  of  radlational  cooling  In  a 
BL.  He  was  motivated  by  the  observation  (Lenschow,  1970)  that  eddy 
flux  decreases  markedly  with  height  throughout  the  BL.  He  found  that 
for  typical  circumstances  In  real  atmosphere,  no  equilibrium  BL  exists 
resulting  from  a balance  between  eddy  heat  flux  convergence  and  di- 
vergence of  radiative  flux,  because  the  latter  is  relatively  too  weak 
and  other  factors  predominate.  However,  he  found  that  if  water  vapor 
is  mostly  concentrated  near  the  surface  (in  fog,  say),  an  equilibrium 
BL  is  possible  due  to  existence  of  a nearly  adiabatic  layer  overlying 
a strongly  superadiabatic  layer.  The  rapid  decrease  of  sensible  heat 
flux  with  height  is  compensated  by  corresponding  radlational  cooling. 
Tennekes  (1973)  suggests  that  since  radiative  heat  flux  tends  to  de- 
crease the  net  heating  rate  of  the  BL,  it  can  be  parameterix&l  by 
making  a small  reduction  in  the  surface  heat  flux. 

GENERAL  REMARKS 

This  chapter  has  discussed  the  determination  of  BL  height,  h,  for 
both  unstable  and  stable  conditions,  which  are  respectively  character- 
ized by  upward  and  downward  surface  heat  flux.  We  have  considered 
both  diagnostic  and  prognostic  methods  for  calculating  h.  We  have 
also  seen  that  diagnostic  methods  are  inadequate  for  describing  an 
evolving  BL.  As  regards  prognostic  equations,  we  have  seen  that  the 
large-scale  vertical  velocity,  and  the  entrainment  at  the  BL,  exercise 
a determining  influence  on  an  unstable  cloud-free  BL.  If  the  BL  has 
cloud  layer  embedded  at  its  top,  the  radiative  heat  flux  also  signifi- 
cantly affects  h.  There  is  some  doubt  regarding  the  use  of  a prognostic 
equation  for  a stable  BL  height. 

The  foregoing  discussions  have  tacitly  assumed  that  the  height  of 
the  thermal  (mixed)  BL  is  equal  to  the  height  of  momentum  as  well  as 
moisture  layers.  Not  many  studies  have  considered  prediction  of  mo- 
mentum BL  height  independently  of  the  thermal  BL  height.  However, 
Deardorff  (1973a)  has  studied  the  effects  of  entrainment  across  a 
velocity  discontinuity  at  the  top  of  a convectively  mixed  BL  having  a 
wind  shear  above.  He  found  that  the  entrainment  produces  large  values 
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of  negative  stress  in  the  upper  two-thirds  of  a BL.  For  a positive 
wind  shear  above  the  BL,  he  obtained  large  positive  values  for  Reynolds 
stress  in  the  upper  portion  of  the  BL. 
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Chapter  S 

ATMOSPHERIC  BOUNDARY  LAYER  AT  LOW  LATITUDES 
GENERAL  FEATURES 

The  dynamics  of  BL  flow  at  low  latitudes  is  not  well  understood 
because  of  the  simultaneous  Importance  of  momentum  transport  by  both 
turbulent  and  large-scale  motions.  Various  observational  and  theo- 
retical studies  have  revealed  the  following  salient  features  of  the 
atmospheric  BL  at  low  latitudes: 

1.  Momentum  advections  are  important  mainly  equatorward  from 
the  Intertropical  Convergence  Zone  (ITCZ)  (Janota,  1971) . 

2.  The  low-level  wind  vector  generally  rotates  with  height,  and 
this  is  mainly  controlled  by  the  height  variation  of  the  horizontal 
pressure  gradient  (Gray , 1968;  Es toque,  1970). 

3.  Advective  accelerations,  associated  with  the  rapid  latitud- 
inal variation  of  coriolis  parameter,  are  important  in  the  BL  between 
the  equator  and  a narrow  latitudinal  transition  (Mahrt,  1972a). 

4.  BL  depth  is  relatively  thin  and  varies  slowly  with  latitude 
even  near  the  equator;  surface  cross-isobar  angle  increases  towards 
the  equator  and  the  cross-isobar  flow  increases  with  height  through- 
out the  BL  (Mahrt,  1972a;  Kraus,  1972). 

5.  Moisture,  in  general,  is  predominantly  confined  to  the  lowest 
levels  near  the  surface,  and  thus  the  BL  is  essential  in  the  thermo- 
dynamics of  (low-latitude)  circulation  systems.  This  was  pointed  out 
by  Charney  and  Eliassen  (1964)  who,  via  the  C1SK  theory,  postulated 
that  the  BL  provides  moisture  convergence  for  latent  heat  release  in 
the  free  atmosphere,  which  in  turn  enables  the  latter  to  induce  upon 
the  boundary  layer  a pressure  field  conducive  to  the  low-level  flow 
convergence.  However,  the  application  of  C7SK  theory  has  suffered 
because  of  lack  of  knowledge  of  BL  dynamics  in  the  low  latitudes,  par- 
ticularly advective  accelerations. 

6.  The  ITCZ,  in  a climatological  sense,  is  a narrow  east-west 
band  of  vigorous  cumulonimbus  convection  and  heavy  precipitation  that 
forms  along  the  equatorward  boundary  of  the  trade-wind  regimes.  (The 
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picture  la  more  complex  whan  viewed  at  ony  particular  instant  in  tine.) 
ITCZ  is  a unique  convective  feature  directly  related  to  the  larger- 
aoale  circulations  in  the  ataosphsre  and  also  to  the  properties  of  the 
oceanic  mixed  layer  (i.e,  the  sea  surface  temperature)  (Pike,  1970). 
Also,  because  of  the  persistnnt  cloudiness  associated  with  ITCZ,  radia- 
tion fluxes  are  important  even  though  the  meridional  extent  of  ITCZ  is 
small. 

7.  A "critical"  latitude  concept  has  been  put  forward  to  explain, 
theoretically,  the  development  of  ITCZ  or  a mean  convergence  sone 
(Holton  et  al. , 1971;  Mahrt,  j.972b).  This  concept  assumes  that  the 
depth  of  the  Bkman  layer  depends  on  the  frequency  of  the  wave  disturb- 
ances in  the  BL,  with  a singularity  at  that  latitude  ct  which  wave 
frequency  is  equal  to  the  coriolis  frequency.  The  range  of  "critical 
latitude"  corresponding  to  this  frequency  band  is  6°  to  7°,  which  is 
dose  to  the  most  frequent  ITCZ  position. 

8.  The  Idea  of  critical  latitude  also  has  been  suggested  by  ob- 
servations, though  on  different  considerations  from  those  noted  above. 
For  example,  Janota  (1968),  on  the  basis  of  detailed  observational 
analysis,  pointed  out  that  flow  characteristics  change  near  4°N  and 
that  the  Ekman  character  of  the  BL  vanishes  south  of  that  latitude. 

9.  There  have  been  various  attempts  to  replace  Ekman-type  param- 
eterization of  BL  convergence  with  the  critical  latitude  concept 
(Yamasaki,  1971;  Hayashi,  1971).  They  all  treat  a BL  with  a top  at 

or  around  900  mb.  However,  observations  (Wallace,  1971;  Reed  and 
Recker,  1971)  have  shown  that  only  10  to  30  percent  of  the  synoptic 
scale  convergence  takes  place  below  900  mb. 

10.  BL  depth  in  the  neighborhood  of  ITCZ  has  no  clear  definition 
because,  with  increasing  convection,  clouds  penetrate  the  stable  in- 
version layer  and  the  fluxes  generated  at  the  surface  may  reach  up  to 
the  high  troposphere.  In  these  regions  (low  latitude),  convective 
activity  generates  a vertical  circulation  that  couples  sea-surface 
processes  with  the  upper  levels  of  the  atmosphere. 

11.  Budget  studies  of  cloud  clusters  embedded  in  trade  winds  and 
ITCZ  (Janota,  1971;  Yanai  et  al.,  1973,  and  others)  have  shown  that 
mean  large-scale  vertical  velocity  is  directed  upward  throughout  the 
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entire  troposphere.  Because  there  is  more  mass  transport  in  clouds 
as  compared  with  the  mean  large-scale  mass  flux,  there  is  a compensat- 
ing sinking  motion  in  clear  areas  between  clouds  that  reaches  its 
maximum  at  cloud  base  (Ogura  and  Cho,  1973;  Yanal  et  al.,  1973).  Thus 
for  strongly  convective  situations  (in  low  latitudes),  the  stable  tran- 
sitional layer  just  below  cloud-base  level  acts  as  an  upper  barrier  for 
turbulent  transports  in  clear  areas  between  clouds. 

PARAMETERIZATION  OF  THE  BOUNDARY  LAYER  IN  LOW-LATITUDE  REGIONS 

It  is  now  well  known  that  low-latitude  regions  are  characterized 
by  deep  convection  of  subsynoptic  scale.  Thus,  if  that  convection 
significantly  influences  large-scale  circulations,  and  if  the  processes 
in  the  BL  have  a determining  influence  on  the  development  or  decay  of 
the  deep  convective  systems,  then  it  is  obvious  that  any  large-scale 
model  must  include  parameterization  of  BL  processes  in  terms  of  large- 
scale  variables. 

Several  different  types  of  BL  have  been  devised — for  instance, 
ir  those  based  on  K-theory  (Es toque  and  Bhumralkar,  1969;  Yamamoto  and 

S' 

Shimanuki,  1966,  etc.),  on  similarity  theory  (Kazanski  and  Monin,  1961; 

i 

Blackadar  and  Tennekes,  1968;  Zilitinkevich,  1969,  1970,  and  others),  ! 

and  so-called  "entity"  models  (Deardorff,  1973;  Tennekes,  1973;  Orlanski  ! 

jr-  v 

et  al.,  1974,  and  others).  But  none  of  these  parameterization  schemes 
‘ f 

is  specifically  designed  to  treat  the  problem  of  the  low-latitude  un-  ] 

I stable,  baroclinic,  Inhomogeneous,  and  transient  atmospheric  BL  com- 

monly found  in  these  regions.  Furthermore,  simple  application  of  the 
existing  models  is  highly  unlikely  to  be  feasible,  at  least  not  until 

thev  are  improved  and  refined  on  the  basis  of  a better  physical  under- 

r- 

| standing  of  BL  processes  in  low  latitudes. 

■ With  this  in  view,  the  Global  Atmospheric  Research  Programme 

(GARP)  has  undertaken  a Boundary-Layer  Subprogramme  (BLSP)  for  the 
GARP  Atlantic  Tropical  Experiment  (GATE),  whose  objectives  have  been 
described  in  GATE  Report  No.  e>,  December  1973.  The  basic  approach  of 
the  BLSP  is  to  use  existing  parameterization  schemes  as  a basis  for 
the  experimental  programme.  One  goal  is  to  improve  those  models  by 
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• Testing  the  assumptions  of  each  scheme,  and 

• Determining  the  importance  of  physical  processes  not 
yet  incorporated  in  the  models. 

It  is  evident  that  any  model  of  the  low-latitude  BL  must  neces- 
sarily consider  the  problem  of  the  dynamics  at  low  latitudes,  where 

the  coriolis  parameter  changes  significantly,  accelerations  affect  | 

the  boundary  structure,  and  cold  ocean  water  generates  strong  baro-  | 

clinity  Uioeber,  1974).  Also,  the  most  important  question  relating  | 

to  the  interaction  between  the  mixed  layer  (usually  found  in  oceanic  § 

* 

low  latitudes)  and  the  cloud  layer  has  to  be  considered.  j 

Recently,  Peardorff  (unpublished)  has  developed  a parameteriza-  I 

I 

tion  scheme  specifically  for  the  low-latitude  BL.  His  scheme  is  based  ] 

cn  the  general  method  suggested  by  Deardorff  (1972r’v,  Betts  (1973),  | 

Tennekes  (1973),  and  Carson  (1973)  for  the  clear-skies  case  and  by  j 

Lilly  (1968)  for  the  case  of  solid  stratocumulus  overcast.  However, 

he  has  extended  this  method  to  apply  to  conditions  of  partial  cumulus 

cloud  cover  (with  bases  at  the  top  of  the  transition  layer),  which  is  I 

-3 

more  usually  prevalent  in  low  latitudes.  He  has  expressed  the  entrain-  = 

I 

ment  velocity  (W  in  Chap.  4),  large-scale  vertical  velocity  at  the  ) 

•2  a. 

top  of  the  BL  (W  ) , and  net  radiative  heat  flux  from  cloud  edges  and  | 

” I 

tope  tn  terms  of  fractional  cloud  cover.  Tests  of  his  modified  param-  * 

eterization  scheme  have  generally  given  realistic  results;  e.g.,:  j 


• Cloud-induced  subsidence  controls  the  growth  of  the  mixed 
BL; 

• Cloud-base  height  stays  just  above  the  top  of  the  mixed 
layer; 

• Sea  surface  temperature  and  large-scale  vertical  motion 
variations  control  rhe  height  of  the  mixed  layer,  cloudi- 
ness, and  other  properties; 

• A vertical  velocity  (subsidence)  of  2 cm/sec  at  a height 
of  1 km  is  sui.icient  to  cause  the  complete  disappearance 
of  convective  clouds. 


Some  of  Deardorff 's  findings.,,  especially  the  one  regarding  the 
role  of  sinking  motion,  have  also  been  obtained  through  much  Simpler 
models  that  treat  the  mixed  layer  as  a slab  (Geisler  and  Kraus,  1969). 
The  "slab"  model  computations  have  shown  that  the  mixed  layer  continues 
to  deepen  unless  the  atmosphere  counteracts  entrainment  by  subsidence. 
The  larger  the  sinking  motion,  the  smaller  the  buoyancy  forces  and  the 
shallower  the  mixed  layer.  Large  subsidence  (typical  of  low- latitude 
regions)  also  results  in  a large  deviation  of  the  BL  wind  from  the 
gcostrophic  wind.  The  model  equations  of  Geisler  and  Kraus  include  a 
drag  coefficient  whose  increased  value  causes  the  mixed  layer  to  deepen 
because  of  the  enhanced  generation  of  turbulent  energy.  This  in  turn 
causes  enhanced  entrainment  of  potentially  warm  air  with  a geostrophic 
value  of  momentum  from  above.  However,  because  of  increased  air-sea 
interaction  (caused  by  increased  drag  coefficient)  the  model,  in  sub- 
stance, yields  larger  ageostrophlc  wind  and  lower  temperatures  in  the 
mixed  layer. 
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Chapter  6 

ATMOSPHERIC  BOUNDARY  LAYER  (ABL)  OVER  OCEANS 


The  underlying  surface  of  the  atmosphere  strongly  affects  atmos- 
pheric processes.  Heat  transfer  and  friction,  both  highly  variable 
in  space  and  time,  are  especially  significant.  The  atmospheric  bound- 
ary layer  (ABL)  is  therefore  important  in  studies  of  general  circula- 
tion, particularly  the  ABL  over  oceans,  which  cover  71  percent  of  the 
earth's  surface  and  predominate  in  tropical  low- latitude  regions. 
Consequently,  Chap.  5 is  by  and  large  applicable  to  marine  boundary 
layers  also;  this  chapter,  however,  deals  exclusively  with  BLs  over 
the  sea. 

The  following  are  the  two  essential  physical  differences  between 
ABL  conditions  over  the  ocean  and  over  land. 

1.  Roughness . For  continents,  whose  land  surface  is  almost 
totally  fixed,  we  can  consider  roughness  as  characteristic 
of  the  surface  itself  and  independent  of  the  flow.  The 
rqughness  of  oceans  is  determined  by  interactions  of  the 
turbulent  motions  of  air  and  water  and  is  therefore  vari- 
able (see  Eqs.  (2.9)  and  (2.10)  in  Chap.  2). 

2.  Storage  of  heat.  Compared  with  temperature  variations  over 
land,  the  temperature  of  water  surface  does  not  change 
appreciably.  Also,  the  dominant  daily  variation  between 
strong  instability  and  stability  over  land  is  nearly  absent 
over  sea;  the  stability  or  instability  of  the  marine  BL  is 
determined  predominantly  through  advection. 

It  is  difficult  to  describe  the  turbulent  structure  in  the  marine 
ABL  because  the  sea  surface  permits  only  a time-space  statistical  de- 
scription, and  atmospheric  turbulence  and  the  state  of  the  sea  interact. 
The  state  of  the  sea  cannot  be  described  simply  as  a local  phenomenon 
due  to  waves  '.raveling  out  of  the  generating  area  into  other  regions 
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far  from  their  origin.  In  fact,  some  scientists  (Hasse,  1970)  believe 
that  the  above  interactions  make  it  hopeless  to  try  to  parameterize 
turbulent  transports  in  a marine  ABL. 

However,  Clark  (1970)  suggested  a simplified  approach  to  param- 
eterize an  oceanic  BL.  He  avoided  the  difficulties  of  modeling  air- 
sea  interaction  by  assuming  a constant  sea-surface  temperature  and  a 
variable  roughness  parameter,  zq.  Otherwise,  his  model  equations  were 
the  same  ones  he  used  for  the  ABL  over  land  (see  Chap.  3 above). 

Honin  and  Zilitinkevich  (1970)  considered  the  parameterization  of  a 
marine  ABL  with  respect  to  large-scale  atmospheric  processes  on  the 
assumption  that  the  wind  and  wave  fields  are  mutually  consistent 

(adapted) . They  defined  an  analog  of  local  Reynolds  number  (R  ) for 

2 W 
the  water  surface  by  using  UA/g  as  a scale  for  measuring  height;  and 

instead  of  considering  the  conventional  roughness  length,  the>  used 

an  effective  roughness  length  given  by 


u: 


CF<V 


where  C„  is  a coefficient  depending  on  R . Similarly,  they  expressed 

r W 

differences  of  temperature  (A8)  and  moisture  (Aq)  between  the  surface 
value  and  the  mean  value  (within  the  BL)  as 


A0  - 9,Cf0(Ru> 

4a  * q*CFq  <Rw> 

Cpg,  being  coefficients  depending  on  Bw<  They  regarded  the  sea- 
surface  water  temperature  as  constant,  on  the  consideration  that  oceans 
have  extremely  large  heat  inertia.  Perhaps  they  did  so  because  they 
could  not  determine  sea-surface  water  temperature  by  solving  a heat 
balance  equation  at  the  sea  surface.  Recognizing  that  in  reality  the 
knowledge  of  coefficients  Cy,  CFQ,  Cpq  is  very  poor  indeed,  they, 
following  Chamock  (1955),  specified 
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and 


CF  - 0.035 


0 


Pandolfo  (1970)  has  described  a nonlinear  numerical  model  of  the 
atmosphere-ocean  BL  that  represents  a complex  local  theory  for  the 
study  of  the  vertical  structure  in  the  marine  ABL.  However,  it  is 
necessary  to  prescribe  horizontal  gradients  if  the  model  is  to  incor- 
porate the  effects  of  horizontal  advection.  The  model  uses  eddy  ex- 
change coefficients  (for  dependent  variables)  that  are  functions  of 
vertical  gradients  of  velocity,  temperature  (both  in  atmosphere  and 
ocean) , humidity  (in  atmosphere) , and  salinity  (in  ocean) . The  model 
equations  consider  the  effect  of  ocean  currents  on  the  momentum  source 
terms. 

Radiative  flux  convergence  is  incorporated  as  a temperature  source 
term  for  computing  both  the  air  and  water  temperatures.  The  most  im- 
portant deficiency  of  the  model  appears  to  be  requirement  of  some  ex- 
plicit form  of  stability  and  wind-wave  dependence  of  eddy  coefficients, 
which  is  not  known  at  all. 

Hasse  (1970)  has  also  parameterized  stress,  heat,  and  moisture 
fluxes  in  the  ABL  at  sea  by  using  a drag  coefficient  dependent  upon 
stability.  However,  this  parameterization  holds  only  as  long  as  there 
is  no  damping  of  turbulence  due  to  strong  stability,  which  can  occur 
over  inland  seas  and  in  coastal  regions  with  cold  water.  Also,  drag 
coefficients  are  uncertain  for  higher  wind  velocities  and  thus  the 
determination  of  transports  in  strong  wind  regimes  is  uncertain. 
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Chapter  7 

HIGHER-MOMENT  APPROACH  TO  BOUNDARY-LAYER  MODELING 


GENERAL  DESCRIPTION 

A three-dimensional  numerical  model  of  Deardorff  (1972b)  has 
shown  considerable  capability  of  supplying  many  fine  details  of  BL 
turbulence  structure.  The  model  parameterizes  the  subgrid  Reynolds 
stresses  by  means  of  nonlinear  eddy  coefficients— a technique  proven 
Inadequate  in  the  treatment  of  a stably  stratified  layer  overlying  a 
well-mixed  convective  layer.  We  have  seen  that  the  magnitude  of  the 
eddy  coefficient  for  heat  (in  these  conditions)  proves  too  large  near 
the  top  of  the  BL,  and  thus  excessively  smoothes  out  the  temperature 
jump  so  evident  in  natural  situations.  This  problem  has  been  handled 
by  Wilhelmson  aid  Ogura  (1972),  who  use  an  eddy  coefficient  term  in 
the  equation  of  temperature  deviation  from  the  mean  state,  but  apply 
no  diffusion  to  the  mean  state.  There  are  some  other  limitations  and 
uncertainties  in  using  the  eddy  coefficient  to  parameterize  subgrid 
scale  fluctuations;  above  the  surface  layer,  for  example,  the  magni- 
tude of  the  eddy  coefficieut  and  its  dependence  upon  stability  are 
generally  not  known.  (Other  limitations  were  discussed  in  Chap.  2.) 

To  overcome  the  problems  associated  with  lowest-order  (K-theory) 
closure  assumptions,  a higher-moment  model  has  been  suggested  by 
Donaldson  (1972).  His  approach  is  to  use  subgrid  transport  equations , 
and  thereby  obviate  the  need  to  parameterize  subgrid  (Reynolds)  stresses 
and  fluxes.  In  this  technique,  mean  motion  and  turbulence  are  separ- 
ated in  a manner  usually  employed  in  treating  BL  flow  but,  in  addition, 
the  equations  for  mean  velocity  and  temperature  are  combined  with 
equations  for  higher-order  statistics,  etch  as  variance  and  covariance, 
and  the  resulting  equations  contain  triple  products.  It  may  be  recog- 
nized that  for  a detailed  3-D  numerical  study  of  the  BL,  the  use  of 
higher-moment  theory  involves  equations  for  15  independent  variables 
(9  components  of  velocity  gradient  tensor,  3 components  of  the  temper- 
ature gradient,  and  3 components  of  the  moisture  gradient).  These  15 
equations  are  in  addition  to  a set  of  three  velocity  component. 
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temperature,  moisture,  and  pressure.  (For  typical  examples  of  equa- 
tions for  subgrid  Reynolds  stresses  or  fluxes,  see  Deardorff  (1973; 
1974),  Donaldson  (1973).) 

Experimentation  with  higher-moment  closure  theory  is  in  its  very 
early  stages,  and  most  of  the  studies  have  used  only  second-moment 
equations  for  generalized  subgrid  stresses  and  fluxes.  Several  prob- 
lems must  be  resolved  before  this  technique  can  be  used  successfully 
in  BL  research — especially  for  the  stable  BL.  Some  of  the  problems 
are: 


• Ad  hoc  relationships  between  third-order  terms  and 
gradients  of  second-order  terms  are  assumed.  As  a 
consequence,  in  some  models  (Donaldson,  1972),  more 
undetermined  quantities  have  to  be  Introduced  with 
the  equations. 

• Closure  assumption  uses  a height-dependent  mixing 
length  that  must  be  assumed  (Donaldson,  1972).  This 
assumption  has  been  avoided  in  some  cases  by  taking 

a representative  grid  scale  (A)  as  the  relevant  length 
scale  in  determining  the  magnitude  of  subgrid  scale 
quantities.  In  doing  so,  however,  the  constants  asso- 
ciated with  terms  stemming  from  pressure  fluctuation 
correlations  must  depend  upon  the  grid  aspect  ratios 
Ax/Az,  Ay/Az,  when  the  latter  are  not  unity  (Deardorff, 
1973).  Presently,  not  much  is  known  about  this 
dependence. 

• The  15  equations  for  stress  and  fluxes  require  computer 
storage  at  each  grid  point  and  for  at  least  one  past 
time  increment.  The  complexity  of  these  equations  in 
finite  difference  form  increases  the  burden  on  com- 
puter requirements.  Deardorff  (1974)  indicates  that 
his  3-D  model  using  higher-moment  closure  assumptions 
consumes  2-1/2  times  as  much  computer  time  as  it  would 
if  eddy  coefficients  were  used  Instead.  To  reduce 
computer  requirements,  a steady-state  nonadvecting 


diifeil 


-76- 


observation  of  the  equations  has  been  suggested*  but 
this  version  is  not  self-consistent  and  tensor- 
invariant  . 

APPLICATION  OF  HIGHER-ORDER  MODELING  TO 
THE  ATMOSPHERIC  BOUNDARY  LAYER 

Current  experimentation  in  applying  the  invariant  higher-order 
(second-order)  closure  model  of  turbulent  shear  flow  to  the  ABL  is 
still  in  its  early  stages.  So  far  Donaldson  (1972,  1973),  Deardorff 
(1973,  1974),  and  Wyngaard  et  al.  (1974)  have  published  their  research 
in  this  connection.  Deardorff  (1974)  has  studied  a heated  (well-mixed) 
BL  numerically  in  a 3-D  model  using  64,000  grid  points  (within  a volume 
5 km  x 5 km  x 2 km)  by  using  subgrid  transport  equations  in  place  of 
eddy  coefficient  formulations.  His  results  compare  very  well  with 
Wangara  data  from  southeast  Australia  (Clark  et  al.,  1971).  This 
model  suggests  that  the  momentum  BL  coincides  with  the  mixed  (thermal) 
BL  during  the  hours  of  solar  heating  of  the  surface;  in  other  words, 
he  finds  invalid  the  concept  that,  within  the  mixed  layer,  the  stress 
vanishes  at  the  lowest  height  at  which  the  corresponding  wind-shear 
component  vanishes.  Wyngaard  et  al.  (1974)  also  use  a higher-order- 
closure  model  to  model  the  convective  BL.  They  have  shown  that  al- 
though coriolis  forces  cause  large  production  rates  of  shear,  the 
mechanism  associated  with  mean  wind  shear  prevents  the  stress  level 
from  becoming  large.  They  aver  that,  in  absence  of  thermal  wind,  the 
stress  problems  are  essentially  linear  regardless  of  wind  direction. 

Higher-order  closure  models  have  not  yet  been  applied  in  atmos- 
pheric general  circulation  models,  essentially  because  of  prohibitive 
computer  requirements. 
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